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1 How to model galactic chemical evolution 

Before going into the detailed chemical evolution history of the Milky Way and 
its satellites, it is necessary to understand how to model, in general, galactic 
chemical evolution. The basic ingredients to build a model of galactic chemical 
evolution can be summarized as : 

• Initial conditions; 

• Stellar birthrate function (the rate at which stars are formed from the gas 
and their mass spectrum); 

• Stellar yields (how elements are produced in stars and restored into the 
interstellar medium); 

• Gas flows (infall, outflow, radial flow). 

When all these ingredients are ready, we need to write a set of equations 
describing the evolution of the gas and its chemical abundances which in- 
clude all of them. These equations will describe the temporal variation of the 
gas content and its abundances by mass (see next sections). The chemical 
abundance of a generic chemical species i is defined as: 

* = —■ (i) 

M 

lvl gas 

According to this definition it holds: 

£*' = 1 ' (2) 

i— l,n 

where n represents the total number of chemical species. Generally, in theoret- 
ical studies of stellar evolution it is common to adopt X, Y and Z as indicative 
of the abundances by mass of hydrogen (H), helium (He) and metals (Z), re- 
spectively. The baryonic universe is madeup mainly of H and some He while 
only a very small fraction resides in metals (all the elements heavier than 
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He), roughly 2%. However, the history of the growth of this small fraction of 
metals is crucial for understanding how stars and galaxies were formed and 
subsequently evolved; and last but not least, because human beings exist only 
because of this small amount of metals! We will focus then our attention is 
studying how the metals were formed and evolved in galaxies, with particular 
attention to our own Galaxy. 

1.1 The initial conditions 

The initial conditions for a model of galactic chemical evolution consist in es- 
tablishing whether :a) the chemical composition of the initial gas is primordial 
or pre-enriched by a pre-galactic stellar generation; b) the studied system is 
a closed box or an open system (infall and/or outflow). 

1.2 Birthrate function 

The birthrate function, can be defined as: 

B(M, t) = ip(t)ip(m) (3) 

where the quantity: 

%l>{t) = SFR (4) 

is called the star formation rate (SFR), namely the rate at which the gas is 
turned into stars, and the quantity: 

<p(m) = IMF (5) 

is the initial mass function (IMF), namely the mass distribution of the stars 
at birth. 

The star formation rate 

The most common parametrization of the SFR is the Schimdt (1959) law: 

m = va k gas , (6) 

where k=l-2 with a preference for k = 1.4 ± 0.15, as suggested by Kennicutt 
(1998a) for spiral disks (see Figure 1), and v is a parameter describing the 
star formation efficiency, in other words, the SFR per unit mass of gas, and 
it has the dimensions of the inverse of a time. Other physical quantities such 
as gas temperature, viscosity and magnetic field are usually ignored. 

Other common parametrizations of the SFR include a dependence on the 
total surface mass density besides the surface gas density: 
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Fig. 1. The SFR as measured by Kennicutt (1998a) in star forming galaxies. The 
continuous line represents the best fit to the data and it can be achieved either with 
the SF law in eq. (6) with k = 1.4 or with the SF law in eq. (9). The short, diagonal 
line shows the effect of changing the scaling radius by a factor of 2. Figure from 
Kennicutt (1998a). 



as suggested by observational results of Dopita & Ryder (1994) and taking 
into account the influence of the potential well in the star formation process 
(i.e. feedback between SN energy input and star formation, see also Talbot & 
Arnett 1975). Other suggestions concern the star formation induced by spiral 
density waves (Wyse & Silk 1998) with expressions like: 

m = vVWRT 1 *™, (8) 

or 

lp(t) = 0.017 Qgasdgas OC R' 1 O" gas (9) 

with fig as being the angular rotation speed of gas (Kennicutt 1998a). Also 
this law provides a good fit to the data of Figure 1. 

The initial mass function 

The most common parametrization of the IMF is a one-slope (Salpeter 1955) 
or multi-slope (Scalo 1986,1998; Kroupa ct al. 1993; Chabricr 2003) power 
law. The most simple example of a one-slope power law is: 



ip(m) = am 



(10) 
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generally denned in a mass range of 0.1-100 M Q , where a is the normalization 
constant derived by imposing that mip(m)dm = 1. 

The Scalo and Kroupa IMFs were derived from stellar counts in the solar 
vicinity and suggest a three-slope function. Unfortunately, the same analysis 
cannot be done in other galaxies and we cannot test if the IMF is the same 
everywhere. Kroupa (2001) suggested that the IMF in stellar clusters is a 
universal one, very similar to the Salpetcr IMF for stars with masses larger 
than O.5M . In particular, this universal IMF is: 

x x = 0.3 for 0.08 < M/M Q < 0.50 



x 2 = 1.3 for M/Mq > 0.5 (11) 

However, Weidner & Kroupa (2005) suggested that the IMF integrated 
over galaxies, which controls the distribution of stellar remnants, the number 
of SNe and the chemical enrichment of a galaxy is generally different from 
the IMF in stellar clusters. This galaxial IMF is given by the integral of the 
stellar IMF over the embedded star cluster mass function which varies from 
galaxy to galaxy. Therefore, we should expect that the chemical enrichment 
histories of different galaxies cannot be reproduced by an unique invariant 
Salpeter-like IMF. In any case, this galaxial IMF is always steeper than the 
universal IMF in the range of massive stars. 



How to derive the IMF 

We define the current mass distribution of local Main Sequence (MS) stars as 
the present day mass function (PDMF), n(m). Let us suppose that we know 
n(m) from observations. Then, the quantity n(m) can be expressed as follows: 
for stars with initial masses in the range 0.1-1.0 M Q which have lifetimes larger 
than a Hubble time we can write: 

rtG 

n(m) = / if(m)tp(t)dt (12) 
Jo 

where tc ~ 14 Gyr (the age of the Universe). The IMF, ip(m), can be taken 
out of the integral if assumed to be constant in time, and the PDMF becomes: 

n(m) = ip(m) < ip > t G (13) 

where < ip > is the average SFR in the past. 

For stars with lifetimes negligible relative to the age of the Universe, 
namely for all the stars with m > 2M Q , we can write: 

n(m) = / tp(m)ip(t)dt, (14) 
Jt G -T m 
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where r m is the lifetime of a star of mass m. Again, if we assume that the 
IMF is constant in time we can write: 

n(m) = ip(m)i>(t G )T m (15) 

having assumed that the SFR did not change during the time interval between 
(to — T m) and to- The quantity ip(tc) is the SFR at the present time. 

We cannot derive the IMF betwen 1 and 2 M Q because none of the previous 
semplifying hypotheses can be applied. Therefore, the IMF in this mass range 
will depend on a quantity, b(t G ): 

«»>-|^ c« 

Scalo (1986) assumed: 

0.5 < b(t G ) < 1.5 (17) 

in order to fit the two branches of the IMF in the solar vicinity. In Figure 
2 we show the differences between a single-slope IMF and multi-slope IMFs, 
which are preferred according to the last studies. 




Fig. 2. Upper panel: different IMFs. Lower panel: normalization of the multi-slope 
IMFs to the Salpeter IMF. Figure from Boissier & Prantzos (1999). 
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1.3 Stellar yields 

The stellar yields, namely the amount of newly formed and pre-existing ele- 
ments ejected by stars of all masses at their death, represent a fundamental 
ingredient to compute galactic chemical evolution. They can be calculated by 
knowing stellar evolution and nucleosynthesis. 

I recall here the various stellar mass ranges and their nucleosynthesis prod- 
ucts. In particular: 

• Brown dwarfs: are stars with masses M < O.IMq which never ignite H. 
They do not enrich the interstellar medium (ISM) in chemical elements 
but only lock up gas. 

• Low and Intermediate mass stars (0.8 < M/M Q < 8.0). Calculations are 
available from Marigo et al. (1996), van den Hoeck & Groenewegen (1997), 
Forestini & Charbonncl (1997), Marigo (2001), Meynet & Maedcr (2002), 
Ventura ct al. (2002), Siess et al. (2002), Karakas & Lattanzio (2007). 
These stars produce mainly 4 He, 12 C, 14 N plus some CNO isotopes and 
s-process (A > 90) elements. In Figure 3 we show an example of integrated 
yields from stars in this mass range. 

• Massive stars (8 < M/M Q < 40). In the mass range 10-40 M , available 
calculations are from Woosley & Weaver (1995, hereafter WW95), Langer 
& Hcnkcl (1995), Thielemann et al. (1996), Nomoto et al. (1997), Limongi 
& Chieffi (2003), Rauscher et al. (2002), Meynet & Maeder (2002), Nomoto 
et al. (2006), among others. These stars end their life as Type II SNe and 
explode by core-collapse; they produce mainly a-elements (O, Ne, Mg, Si, 
S, Ca), some Fe-peak elements, s-process elements (A < 90) and r-process 
elements. Stars more massive than 40M© can end up as Type Ib/c SNe: 
they are also core-collapse SNe and are linked to 7-ray bursts (GRB). 

• Type la SNe (white dwarfs in binary systems, see later). Calculations are 
available from Nomoto et al. (1997), Iwamoto et al. (1999). They produce 
mainly Fe-peak elements. 

• Very massive objects (M > 100M Q ). Calculations are available from e.g. 
Portinari ct al. (1998), Umcda & Nomoto (2001). They should produce 
mainly oxygen although many uncertainties are still present. 

All the elements with mass number A from 12 to 60 have been formed 
in stars during the quiescent burnings. Stars transform H into He and then 
He into heaviers until the Fe-peak elements, where the binding energy per 
nucleon reaches a maximum and the nuclear fusion reactions stop. H is trans- 
formed into He through the proton-proton chain or the CNO-cycle, then 4 He 
is transformed into 12 C through the triple- a reaction. 

Elements heavier than 12 C are then produced by synthesis of a-particles: 
they are called a-elements (O, Ne, Mg, Si and others). 

The last main burning in stars is the 28 Si -burning which produces 56 Ni, 
which then decays into 56 Co and 56 Fe. Si-burning can be quiescent or explo- 
sive (depending on the temperature). 



Chemical evolution of the Milky Way and its Satellites 7 

Explosive nucleosynthesis occurring during SN explosions mainly produces 
Fe-peak elements. Elements originating from s- and r-processes (with A> 60 
up to Th and U) are formed by means of slow or rapid (relative to the (3- 
decay) neutron capture by Fe seed nuclei; s-processing occurs during quiescent 
He-burning, whereas r-processing occurs during SN explosions. 

In Figures 4, 5, 6, 7 and 8 we show a comparison between stellar yields 
for massive stars computed for different initial stellar metallicities and with 
different assumptions concerning the mass loss. In particular, some yields are 
obtained by assuming mass loss by stellar winds with a strong dependence on 
metallicity (e.g. Maeder, 1992), whereas others (e.g. WW95) are computed 
by means of conservative models without mass loss. One important difference 
arises for oxygen in massive stars for solar metallicity and mass loss: in this 
case, the O yield is strongly depressed as a consequence of mass loss. In fact, 
the stars with masses > 25M Q and solar metallicity lose a large amount of 
matter rich of He and C, thus subctracting these elements to further process- 
ing which would lead to O and heavier elements. So the net effect of mass 
loss is to increase the production of He and C and to depress that of oxy- 
gen (see Figure 9). More recently, Meynet & Mader (2002, 2003, 2005) have 
computed a grid of models for stars with M > 20M Q including rotation and 
metallicity dependent mass loss. The effect of metallicity dependent mass loss 
in decreasing the O production in massive stars was confirmed, although they 
employed significantly lower mass loss rates compared with Maeder (1992). 
With these models they were able to reproduce the frequency of WR stars and 
the observed WN/WC ratio, as was the case for the previous Maeder results. 
Therefore, it appears that the earlier mass loss rates made-up for the omis- 
sion of rotation in the stellar models. On the other hand, the dependence upon 
metallicities of the yields computed with conservative stellar models, such as 
those of WW95, is not very strong except perhaps for the yields computed 
with zero intial stellar metallicity (Pop III stars). 

In Figures 7 and 8 we show the most recent results of Nomoto et al. (2006) 
for conservative stellar models of massive stars at different metallicities. While 
the O yields are not much dependent upon the initial stellar metallicity, as in 
WW95 , the Fe yields seem to change dramatically with the stellar metallicity. 

Type la SN progenitors 

There is a general consensus about the fact that SNela originate from C- 
deflagration in C-0 white dwarfs (WD) in binary systems, but several evo- 
lutionary paths can lead to such an event. The C-dcflagration produces 
~ 0.6 — 0.7M Q of Fe plus traces of other elements from C to Si, as observed 
in the spectra of Type la SNe. 

Two main evolutionary scenarios for the progenitors of Type la SNe have 
been proposed: 

• Single Degenerate (SD) scenario (see Figure 10): the classical scenario of 
Whclan and Iben (1973), recently revised by Han & Podsiadlowsky (2004), 
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Fig. 3. The yields integrated over the Salpeter (1955) IMF of He, C and N produced 
by low and intermediate mass stars as functions of the initial stellar metallicity. 
Different results are compared here: those of RV81 (Renzini & Voli 1981), those of 
HG97 (van den Hoeck & Groenewegen 1997) and those of M2K (Marigo 2001). The 
mixing length parameters (a) adopted by the authors are indicated. Figure from 
Marigo (2001). 



namely C-dcflagration in a C-0 WD reaching the Chandrasekhar mass 
Mch ~ 1-44M after accreting material from a red giant companion. One 
of the limitations of this scenario is that the accretion rate should be de- 
fined in a quite narrow range of values. To avoid this problem, Kobayashi et 
al. (1998) had proposed a similar scenario, based on the model of Hachisu 
et al. (1996), where the companion can be either a red giant or a main 
sequence star, including a metallicity effect which suggests that no Type 
la systems can form for [Fe/H]< —1.0 dex. This is due to the development 
of a strong radiative wind from the C-0 WD which stabilizes the accre- 
tion from the companion, allowing for larger mass accretion rates than the 
previous scenario. The clock to the explosion is given by the lifetime of the 
secondary star in the binary system, where the WD is the primary (the 
originally more massive one). Therefore, the largest mass for a secondary 
is 8M Q , which is the maximum mass for the formation of a C-0 WD. As a 
consequence, the minimum timescale for the occurrence of Type la SNe is 
~ 30 Myr (i.e. the lifetime of a 8M ) after the beginning of star formation. 
Recent observations in radio-galaxies by Mannucci et al. (2005; 2006) seem 
to confirm the existence of such prompt Type la SNe. 
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Fig. 4. The yields of oxygen for massive stars as computed by several authors, as 
indicated in the Figure. None of these calculations takes into account mass loss by 
stellar wind. 
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Fig. 5. The same as Fig. 4 for magnesium. 
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Fig. 7. The O yields as calculated by Nomoto et al. (2006) for different metallicities. 
These calculations do not take into account mass loss by stellar wind. 
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Fig. 8. The same as Figure 7 for Fe. 
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Fig. 9. The effect of metallicity dependent mass loss on the oxygen yield. The 
comparison is between the conservative yields of WW95 for Z=0.001 and Z=0.02 
and the yields with mass loss of Maeder (1992) for the same metallicity. As one can 
see the effect of mass loss for a solar metallicity is a quite important one. 
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The minimum mass for the secondary is O.8M , which is the star with 
lifetime equal to the age of the universe. Stars with masses below this limit 
are obviously not considered. In summary, the mass range for both primary 
and secondary stars is, in principle, between 0.8 and 8M , although two 
stars of O.8M0 are too small to give rise to a WD with a Chandrasekhar 
mass, and therefore the mass of the primary star should be assumed to 
be high enough to ensure that, even after accretion from a 0.8M© star 
secondary, it will reach the Chandrasekhar mass. 
• Double Degenerate (DD) scenario: the merging of two C-0 white dwarfs, 
due to loss of angular momentum caused by gravitational wave radiation, 
which explode by C-dcflagration when Mqh is reached (Iben and Tutukov 
1984). In this scenario, the two C-0 WDs should be of ~ 0.7M Q in order 
to give rise to a Chandrasekhar mass after they merge, therefore their 
progenitors should be in the range (5-8)Mq. The clock to the explosion 
here is given by the lifetime of the secondary star plus the gravitational 
time delay which depends on the original separation of the two WDs. The 
minimum timescale for the appearance of the first Type la SNe in this 
scenario is a few million years more than in the SD scenario (e.g. ~ 40 
Myr in Tornambe & Matteucci 1986). At the same time, the maximum 
gravitational time delay can be as long as more than a Hubble time. For 
more recent results on the DD scenario see Greggio (2005). 

Within any scenario the explosion can occur cither when the C-0 WD reaches 
the Chandrasekhar mass and carbon deflagrates at the center or when a mas- 
sive enough helium layer is accumulated on top of the C-0 WD. In this last 
case there is He-detonation which induces an off-center carbon deflagration be- 
fore the Chandrasekhar mass is reached (sub-chandra exploders, e.g. Woosley 
& Weaver 1994). 

While the chandra-exploders are supposed to produce the same nucle- 
osynthesis (C-dcflagration of a Chandrasekhar mass), they predict a different 
evolution of the Type la SN rate and different typical timescales for the SNe 
la enrichment. A way of defining the typical Type la SN timescale is to as- 
sume it as the time when the maximum in the Type la SN rate is reached 
(Matteucci & Recchi, 2001). This timescale varies according to the chosen pro- 
genitor model and to the assumed star formation history, which varies from 
galaxy to galaxy. For the solar vicinity, this timescale is at least 1 Gyr, if the 
SD scenario is assumed, whereas for elliptical galaxies, where the stars formed 
much more quickly, this timescale is only 0.5 Gyr (Matteucci & Greggio, 1986; 
Matteucci & Recchi 2001). 

1.4 Gas flows 

Various parametrizations have been suggested for gas flows and the most 
common is an exponential law for the gas infall rate: 

IRcte- t/T (18) 
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The progenitor of a Type la supernova 
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Fig. 10. The progenitor of a Type la SN in the context of the single-degenerate 
model (Illustration credit: NASA, ESA, and A. Field (STSci)). 



with the timescale r being a free parameter, whereas for the galactic outflows 
the wind rate is generally assumed to be proportional to the SFR: 

WR = -XSFR (19) 

where A is again a free parameter. Both r and A should be fixed by reproducing 
the majority of observational constraints. 
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2 Basic Equations for chemical evolution 
2.1 Yields per Stellar Generation 

Under the assumption of Instantaneous Recycling Approximation (IRA) 
which states that all stars more massive than 1M@ die immediately, whereas 
all stars with masses lower than 1M live forever, one can define the yield per 
stellar generation (Tinsley, 1980): 

1 /- 00 

Vi = 1 _ R J mp im ip(m)dm (20) 

where pi m is the stellar new yield of the element i, namely the newly 
formed and ejected element i by a star of mass m, and ip(m) is the IMF. 
The quantity R is the so-called Returned Fraction: 

(m — M rem )(p(m)dm (21) 

and is the mass fraction of gas restored into the ISM by an entire stellar 
generation. The term fraction derives from the fact that in its definition R is 
divided by J °° mip(m)dm = 1, which is the normalization condition for the 
IMF. Therefore, (1-R) is the fraction of mass locked up in very low mass stars 
and remnants. 



2.2 Analytical models 

Simple Model 

The Simple Model for the chemical evolution of the solar neighbourhood is 
the simplest approach to model chemical evolution. The solar neighbourhood 
is assumed to be a cylinder of 1 Kpc radius centered around the Sun. 
The basic assumptions of the Simple Model are: 

- the system is one-zone and closed, no inflows or outflows are considered, 

- the initial gas is primordial (no metals), 

- IRA holds, 

- the IMF, <p(m), is assumed to be constant in time, 

- the gas is well mixed at any time (instantaneous mixing approximation, 
IMA). 

The Simple Model fails in describing the evolution of the Milky Way 
(G-dwarf metallicity distribution, elements produced on long timescales and 
abundance ratios) and the reason is that at least two of the above assumptions 
are manifestly wrong, epecially if one intends to model the evolution of the 
abundance of elements produced on long timescales, such as Fe. In particular, 
the assumptions of the closed box and the IRA. 
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However, it is interesting to know the solution of the Simple Model and 
its implications. Let Z be the abundance by mass of metals, then we obtain 
the analytical solution of the Simple Model by ignoring the stellar lifetimes: 

Z = yzln(^) (22) 

where G — M gas /M t ot is the gas mass fraction of the system and yz is the 
yield per stellar generation, as defined above, otherwise called effective yield. 
In particular, the effective yield is defined as: 

yz -" = mTJg) (23) 

namely the yield that the system would have if behaving as the simple closed- 
box model. This means that if yz eff > yz, then the actual system has attained 
a higher metallicity at a given gas fraction G. Generally, given two chemical 
elements i and j, the solution of the Simple Model for primary elements (eq.22) 
implies: 

Ji = y± (24) 

which means that the ratio of two element abundances is always equal to 
the ratio of their yields. This is no more true when IRA is relaxed. In fact, 
relaxing IRA is necessary to study in detail the evolution of the abundances 
of single elements produced on long timescales (e.g. Fe, N). 



Analytical models in the presence of outflows 

One can obtain analytical solutions also in the presence of infall and/or outflow 
but the necessary condition is to assume IRA, as well as precise forms for the 
infall and outflow rates. 

Matteucci & Chiosi (1983) found solutions for models with outflow and 
infall and Matteucci (2001) found it for a model with infall and outflow acting 
at the same time. The main assumption in the model with outflow but no infall 
is that the outflow rate is: 



W(t) = -A(l - R)%l>(t) (25) 

where A > is the wind parameter. 
The solution of this model is: 

z= (ITA) /n[(1 + A)G_1_A] (26) 

for A = the equation becomes the one of the Simple Model (eq. 22). 
As one can see from eq. (26), the presence of an outflow decreases the 
effective yield, in the sense that the true yield of a system is lower than the 
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effective yield. Models with galactic winds or outflows in general are suit- 
able for ellipticals, irregulars and for the Galactic halo. A popular analytical 
model with outflow is that suggested by Hartwick (1976) for the evolution of 
the Galactic halo, under the assumption that during the halo collapse stars 
were forming while the gas was dissipating energy and falling into the bulge 
and disk, thus producing a net gas loss from the halo. This hypothesis was 
suggested by the fact that the stellar metallicity distribution of the halo can 
be reproduced only with an effective yield lower than that of the disk. In 
Hartwick's model the ouflow rate is assumed to be simply proportional to the 



which is similar to eq. (25). Hartwick used this model to reproduce the 
metallicity distribution of halo stars and also to alleviate the G-dwarf problem 
in the disk, namely the fact that the Simple Model of chemical evolution 
predicts too many disk stars than observed (see Tinsley 1980 for a review on 
the subject). However, the gas lost from the halo cannot have contributed to 
form the whole disk since the distribution of the specific angular momentum 
of halo and disk stars are quite different, thus indicating that only a negligible 
amount of halo gas can have formed the disk. On the other hand, the similarity 
of the distributions for the halo and bulge indicates that the bulge must have 
formed out of gas lost from the halo (see Wyse & Gilmore 1992). The G- 
dwarf problem is instead easily solved if one assumes that the Galactic disk 
has formed by means of slow infall of extragalactic material, as we will see in 
the next sections. Recently, Hartwick's model has been revisited by Prantzos 
(2003) to interpret the most recent metallicity distribution of halo stars, which 
is quite different with respect to the G-dwarf metallicity distribution in the 
local disk. In particular, the halo metallicity distribution is peaked at around 
[Fe/H]=-1.6 dex, whereas the G-dwarf distribution is peaked at around ~ 
—0.2 dex. Prantzos (2003) suggested that an outflow with A = 8 as well as a 
formation of the halo by early infall are necessary to reproduce the observed 
halo metallicity distribution. In Figure 11 we show the results of Prantzos 
(2003) compared with observations. 

Analytical models in presence of infall 

The solution of the equation of metals for a model without a wind but with 
a primordial infalling material (Za = 0) at a rate: 



SFR: 



W(t) 



(27) 



A(t) = A(l 



mm 



(28) 



and A ^ 1 is : 




(29) 
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Fig. 11. Metallicity distribution for the halo stars. Upper panel: observed and 
predicted metallicity distributions. The models are: pure outflow with IRA (dashed 
curve) , pure outflow without IRA (thin solid curve) and early infall +outflow without 
IRA (thick solid curve). The distribution is on a linear scale. Middle panel: the 
same as above but the distribution is on a logarithmic scale. Lower panel: predicted 
cumulative distributions. Figure from Prantzos (2003) . 

For A — 1 one obtains the well known case of extreme infall studied by Larson 
(1972) whose solution is: 

Z = y z [\-e-^ 1 ^] (30) 

This extreme infall solution shows that when G — » then Z — > yz- The infall 
can solve the G-dwarf problem for disk stars except for the extreme infall 
solution which predicts too few low metallicity stars below [Fe/H] =-1.0 (see 
Tinsley 1980). Finally, we can conclude that the infall is very important for 
explaining both the halo and the disk formation. 
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Analytical models in presence of infall and outflow 

Matteucci (2001) presented an analytical solution for infall and outflow 
present at the same time. The solution refers to the outflow and infall rates 
of eq. (25) and eq. (28), respectively. 
In particular: 

z =yz{i-[(A-\) 

for a primordial infalling gas (Za = 
A > and A ± 1. 

2.3 Detailed numerical models 

Detailed models of galactic chemical evolution require consideration of the 
stellar lifetimes, namely they should relax IRA. However, the majority of them 
still retain the instantaneous mixing approximation (IMA), which assumes 
that the material ejected by stars at their death is instantaneously mixed 
with the surrounding interstellar medium (ISM). This approximation seems 
to be good in the majority of the cases with perhaps the exception of the very 
early phases of galactic evolution. 

The basic equations of chemical evolution follow the evolution of the abun- 
dances of single chemical species and the g whole. 

If (Tj is the surface mass density of an element i, with a gas = ^2 i=1 n o% be- 
ing the total surface gas density and n the total number of chemical elements, 
we can write: 



-(vl-A-l)^- 1 ]*^}, (31) 
= 0). This solution is velid for A > and 



&i(t) = -i>(t)Xi(t) 

rM Bm 
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J M L 

+A / (j)(m) 

JM Bm 
/■0.5 

•[ / f(l)lP(t- T m 2)Qmi{t - T m 2)dj}dm 
i-Mbm 

+B tp(t - T m )Q mi (t - T m )tp(m)dm 

JM Bm 
,-Mu 

+ ip(t - T m )Q mi (t - T m )ip(m)dm 

JM BM 

+X Ai A(t) - Xi(t)W(t) (32) 

for any ggiven chemical element. These equations can be solved only nu- 
merically. The quantities Xi{t) are the abundances as defined in eq. (1). The 
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quantity Q m i contains all the information about stellar evolution and nucle- 
osynthesis: in practice it gives the mass of gas produced and ejected in the 
form of an element i by a star of initial mass m, together with the mass of that 
element which was already present in the star at birth. The various integrals 
represent the rates at which the mass of a given element is restored into the 
ISM by stars of different masses which can evolve into WDs or supernovae 
(II, la, lb). The integral representing the rate of matter restoration by Type 
la SNe is the second one on the right hand side. The quantity A is a constant: 
it is the fraction, in the IMF, of binary systems with those specific features 
required to give rise to Type la SNe, whereas B=l-A is the fraction of all the 
single stars and binary systems in the same mass range of definition of the 
progenitors of Type la SNe (third integral) . The parameter A is obtained by 
imposing that the predicted Type la SN rate reproduces the observed rate at 
the present time (14 Gyr). Values of A =0.05-0.09 are found for the evolution 
of the solar vicinity when an IMF of Scab (1986, 1989) or Kroupa et al. (1993) 
is adopted. If one adopts a flatter IMF such as the Salpeter (1955) one, then A 
is different. The integral of the Type la SN contribution is made over a range 
of mass going from Msm — 3M Q to Mbm — 16A/ , which represents the 
total masses of binary systems able to produce Type la SNe in the framework 
of the single degenerate scenario. There is also an integration over the mass 
distribution of binary systems; in particular, one considers the function f(j) 
where 7 = M ^ M2 , with Mi and M2 being the primary and secondary mass 
of the binary system, respectively (for more details see Matteucci & Greggio 
1986 and Matteucci 2001). The third and fourth integrals represent the rates 
of Type II and Type Ib/c SNe, respectively. The occurrence of Type lb SNe 
seems to be partly related to Wolf-Rayet stars which have original masses 
larger tham 25M© and depends on the mass loss rate which is more active at 
high metallicities. However, it has been proposed that Type lb SNe can also 
originate from massive stars in binary systems. Finally, the functions A(t) and 
W(t) are the infall and wind rate, respectively. 

3 The Milky Way 

We will first analyze the chemical evolution of our Galaxy, the Milky Way. 

3.1 The formation of the Milky Way 
Observational evidence 

The Milky Way galaxy has four main stellar populations: 1) the halo stars 
with low metallicities □ and eccentric orbits, 2) the bulge population with a 



3 the most common stellar metallicity indicator in stars is [Fe/H]= log(Fe/ H), — 
log(Fe/H) e 
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Fig. 12. Schematic edge-on view of the major components of the Milky Way. Illus- 
tration credit from R. Buser, www. astro. unibas.ch/forschung/rb/structure.shtml. 

large range of metallicities and dominated by random motions, 3) the thin disk 
stars with an average metallicity < [Fe/H] >=-0.5 dex and circular orbits, 
and finally 4) the thick disk stars which possess chemical and kinematical 
properties intermediate between those of the halo and those of the thin disk. 
The halo stars have average metallicities of < [Fe/H] >=-1.5 dex and a 
maximum metallicity of ~ —1.0 dex, although stars with [Fe/H] as high as 
-0.6 dex and halo kinematics are observed. The average metallicity of thin 
disk stars is ~ —0.6 dex, whereas the one of bulge stars is ~ —0.2 dex. 

The kinematical and chemical properties of the different Galactic stellar 
populations can be interpreted in terms of the Galaxy formation mechanism. 
Eggen, Lynden-Bell & Sandage (1962), in a cornerstone paper suggested a 
rapid monolithic collapse for the formation of the Galaxy lasting ~ 2 • 10 8 
years. This suggestion was based on a kinematical and chemical study of solar 
neighbourhood stars and the value of the suggested timescale was chosen to 
allow for the orbital eccentricities to vary in a potential not yet in equilibrium 
but sufficiently long so that massive stars forming in the collapsing gas could 
have time to die and enrich the gas with heavy elements. 

Later on, Searle & Zinn (1978) measured Fe abundances and horizontal 
branch morphologies of 50 globular clusters and studied their properties as a 
function of the galactocentric distance. As a result of this, they proposed a 
central collapse like the one envisaged by Eggen et al. , but also that the outer 
halo formed by merging of large fragments taking place over a considerable 
timescale > 1 Gyr. The Searle & Zinn scenario is close to what is predicted by 
modern cosmological theories of galaxy formation. In particular, in the frame- 
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work of the hierarchical galaxy formation scenario, galaxies form by accretion 
of smaller building blocks (e.g. White & Rees, 1978, Navarro & al. 1997). Ob- 
vious candidates for these building blocks are either dwarf spheroidal (dSph) 
or dwarf irregular (dlrr) galaxies. However, as we will see in detail later, the 
chemical composition and in particular the chemical abundance patterns in 
dSphs or dlrrs arc not compatible with the same abundance patterns in the 
Milky Way (see Geisler et al. 2007), thus arguing against the identification 
of the building blocks with these galaxies. On the other hand, very recently, 
Carollo & al. (2007) have obtained medium resolution spectroscopy of 20,336 
stars from the Sloan Digital Sky Survey (SDSS). They showed that the Galac- 
tic halo is divisible into two broadly overlapping structural components. In 
particular, they find that the inner halo is dominated by stars with very ec- 
centric orbits, exhibits a peak at [Fe/H]=-1.6 dex and has a flattened density 
distribution with a modest net prograde rotation. The outer halo includes 
stars with a wide range of eccentricities, exhibits a peak at [Fe/H]=-2.2 dex 
and a spherical density distribution with highly statistically significant net ret- 
rograde rotation. They conclude that most of the Galactic halo should have 
formed by accrection of multiple distinct sub-systems. However, an analysis 
of the abundance ratios of these stars is still missing. 

Theoretical models 

From an historical point of view, the modelization of the Galactic chemical 
evolution has passed through different phases that I summarize in the follow- 
ing: 

• SERIAL FORMATION 

The Galaxy is modeled by means of one accretion episode lasting for the 
entire Galactic lifetime, where halo, thick and thin disk form in sequence 
as a continuous process. The obvious limit of this approach is that it does 
not allow us to predict the observed overlapping in mctallicity between 
halo and thick disk stars and between thick and thin disk stars, but it 
gives a fair representation of our Galaxy (e.g. Matteucci & Francois 1989). 

• PARALLEL FORMATION 

In this formulation, the various Galactic components start at the same time 
and from the same gas but evolve at different rates (e.g. Pardi et al. 1995). 
It predicts overlapping of stars belonging to the different components but 
implies that the thick disk formed out of gas shed by the halo and that 
the thin disk formed out of gas shed by the thick disk, and this is at 
variance with the distribution of the stellar angular momentum per unit 
mass (Wyse & Gilmore 1992), which indicates that the disk did not form 
out of gas shed by the halo. 

• TWO-INFALL FORMATION 

In this scenario, halo and disk formed out of two separate infall episodes 
(overlapping in mctallicity is also predicted) (e.g. Chiappini et al. 1997; 
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Chang et al. 1999; Alibes et al. 2001). The first infall episode lasted no 
more than 1-2 Gyr whereas the second, where the thin disk formed, lasted 
much longer with a timescale for the formation of the solar vicinity of 6-8 
Gyr (Chiappini et al. 1997; Boissierfe Prantzos 1999). 
• STOCHASTIC APPROACH 

Here the hypothesis is that in the early halo phases ([Fe/H] < —3.0 dex), 
mixing was not efficient and, as a consequence, one should observe, in 
low metallicity halo stars, the effects of pollution from single SNe (e.g. 
Tsujimoto et al. 1999; Argast et al. 2000; Oey 2000). These models predict 
a large spread for [Fe/H] < —3.0 dcx for all the a-elements, which is not 
observed, as shown by stellar data with metallicities down to -4.0 dex 
(Cayrel et al. 2004). However, inhomogencities could explain the observed 
spread of s- and r-elements at low metallicities (see later). 

3.2 The two-infall model 

The two-infall model of Chiappini, Matteucci & Gratton (1997) predicts two 
main episodes of gas accretion: during the first one, the halo the bulge and 
most of the thick disk formed, while the second gave rise to the thin disk. In 
Figure 13 we show an artistic representation of the formation of the Milky Way 
in the two-infall scenario. In the upper panel we see the sequence of the forma- 
tion of the stellar halo, in particular the inner halo, following a monolithic-like 
collapse of gas (first infall episode) but with a longer timescale than originally 
suggested by Eggen et al. (1962): here the time scale is 0.8-1.0 Gyr. During the 
halo formation also the bulge is formed on a very short timescale in the range 
0.1-0.5 Gyr. During this phase also the thick disk assembles or at least part of 
it, since part of the thick disk, like the outer halo, could have been accreted. 
The second panel from left to right shows the beginning of the thin disk forma- 
tion, namely the assembly of the innermost disk regions just around the bulge; 
this is due to the second infall episode. The thin-disk assembles inside-out, in 
the sense that the outermost regions take a much longer time to form. This 
is shown in the third panel. In Fig. 13 each panel is connected to temporal 
phases where the Type II and then the Type la SN rates are present. So, it 
is clear that the early phases of the halo and bulge formation are dominated 
by Type II SNe (and also by Type Ib/c SNe) producing mostly a-elements 
such as O and Mg and part of Fe. On the other hand, Type la SNe start to 
be non negligible only after lGyr and they pollute the gas during the thick 
and thin disk phases. The minimum shown in the Type II SN rate is due to a 
gap in the star formation rate occurring as a consequence of the adoption of 
a threshold density in the star formation process, as we will see next (Figure 
14). 

3.3 Detailed recipes for the two-infall model 

The main assumption of this model are: 
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Fig. 13. Artistic view of the two-infall model by Chiappini et al. (1997). The pre- 
dicted SN II and la rates per century are also sketched, together with the fact that 
Type II SNe produce mostly a-elements (e.g. O, Mg), whereas Type la SNe pro- 
duce mostly Fe. (Illustration credit: C. Chiappini, Sky & Telescope, 2004, Vol. 108, 
number 4, p. 32 ). 

• The IMF is that of Scalo (1986) normalized over a mass range of 0.1- 
100M Q . 

• The infall law is: 

A(r, t) = a{r)e- t/TH{r) + &(r)e- ( *-'""» )/tdW (33) 

where A(r,t) = )infaU is the ra te at which the total surface mass 

density changes because of the infalling gas. The quantities a(r) and b(r) 
are two parameters fixed by reproducing the total present time surface 
mass density in the solar vicinity (cr tot = 51 ± 6 Mq pc~ 2 , see Boissier 
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& Prantzos 1999), t max = lGyr is the time for the maximum infall on 
the thin disk, th — 0.8 Gyr is the time scale for the formation of the 
halo thick-disk (which means a total duration of 2 Gyr for the complete 
halo-thick disk formation) and td(t) is the timescale for the formation of 
the thin disk and it is a function of the galactocentric distance (formation 
inside-out, Matteucci & Francois 1989; Chiappini et al. 2001). 
In particular, it is assumed that: 

t d = lM3r(Kpc) - 1.267 (Gyr) (34) 

where r is the galocentric distance. 

• The SFR is the Kcnnicutt law with a dependence on the surface gas density 
and also on the total surface mass density (see Dopita & Ryder 1994). In 
particular, the SFR is based on the law originally suggested by Talbot & 
Arnett (1975) and then adopted by Chiosi (1980): 

^(^) = ^( a(r ^ M) ) (fc " 1) ^(^) fc - (35) 

where the constant v is the efficiency of the SF process, as defined in eq. 
(6), and is expressed in Gyr' 1 : in particular, v — 2Gyr~ 1 for the the halo 
and lGyr^ 1 for the disk (t > lGyr). The total surface mass density is 
represented by a(r, t), whereas a(r Q , t) is the total surface mass density at 
the solar position, assumed to be r Q = 8 Kpc (Reid 1993). The quantity 
<Jgas (r, t) represents the surface gas density. The exponent of the surface 
gas density, k, is set equal to 1.5, similar to what suggested by Kennicutt 
(1998a). These choices for the parameters allow the model to fit very well 
the observational constraints, in particular in the solar vicinity. We recall 
that below a critical threshold for the surface gas density (7M & pc~ 2 for 
the thin disk and 4M©pc~ 2 for the halo phase) we assume that the star 
formation is halted. The existence of a threshold for the star formation has 
been suggested by Kennicutt (1998a, b) and Martin & Kennicutt (2001). 
The predicted behaviour of the SFR, obtained by adopting eq.(35) with 
the threshold is shown in Figure 14. 

• The assumed Type la SN model is the single-degenerate one with the 
recipe first adopted in Greggio & Renzini (1983a) and Matteucci & Greggio 
(1986) and more recently in Matteucci & Recchi (2001). The minimum 
time for the explosion is 30 Myr, whereas the the timescale for restoring 
the bulk of Fe is 1 Gyr, for the SFR adopted in the solar vicinity. It is 
worth recalling that this timescale is not universal since it depends on the 
assumed SNIa progenitor model but also on the assumed star formation 
history. The SN rates in the solar vicinity are shown in Figure 15. 

3.4 The chemical enrichment history of the solar vicinity 

We study first the solar vicinity, namely the local ring at 8 Kpc from the 
galactic center. By integrating eq.(32) without the wind term we obtain the 
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Fig. 14. The SFR in the solar vicinity as predicted by the two-infall model. Figure 
from Chiappini et al. (1997). The oscillating behaviour in the SFR at late times 
is due to the assumed threshold density for SF. The threshold gas density is also 
responsible for the gap in the SFR seen at around 1 Gyr. 

evolution of the abundances of several chemical species (H, D, He, Li, C, N, 
O, a-elements, Fe, Fe-peak elements, s-and r- process elements). In Figure 16 
we show the smallest mass dying at any cosmic time corresponding to a given 
predicted abundance of [Fe/H] in the ISM. This is because there is an age- 
metallicity relation and the [Fe/H] abundance increases with time. We recall 
that, for a generic chemical element i, with abundance Xi, one defines: 

[Xi/H] = log(Xi/H)* - log(X l /H) G , (36) 

where logiXij H)q refers to the solar abundance of the element i. 

The observational constraints 

A good model of chemical evolution should be able to reproduce a minimum 
number of observational constraints and the number of observational con- 
straints should be larger than the number of free parameters which are: th, 
to, ki, k2, v and A (the fraction of binary systems which can give rise to 
Type la SNe). 

The main observational constraints in the solar vicinity that a good model 
should reproduce (see Chiappini et al. 2001, Boissier & Prantzos, 1999 and 
references therein) are: 
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Fig. 15. The Type II and la rate in the solar vicinity as predicted by the two-infall 
model. Figure from Chiappini et al. (1997). The oscillating behaviour of the Type II 
SN rate at late times is due to the assumed threshold density for SF. The threshold 
gas density is also responsible for the gap in the SFR seen at around 1 Gyr. 



• The present time surface gas density: <j gas = 13 ± 3M & pc~ 2 

• The present time surface star density c* = 43 ± 5M©pc~ 2 

• The present time total surface mass density: atot = 51 ± 6Mqpc~ 2 

• The present time SFR: Vo = 2 - 5M Q pc~ 2 Gyr _1 

• The present time infall rate: 0.3 — 1 .5Mqpc~ 2 'Gyr -1 

• The present day mass function (PDMF) 

• The solar abundances, namely the chemical abundances of the ISM at the 
time of birth of the solar system 4.5 Gyr ago as well as at the present time 
abundances 

• The observed [Xi/Fe] vs. [Fe/H] relations 

• The G-dwarf metallicity distribution 

• The age-met allicity relation 

And finally, a good model of chemical evolution of the Milky Way should re- 
produce the distributions of abundances, gas and star formation rate along the 
disk as well as the average observed SNII and la rates (SNII=1.2±0.8 10Ch/r -1 
and SNIa=0.3 ± 0.2 lOOyr-- 1 ). 

The time-delay model 

What we call time-delay model is the interpretation of the behaviour of abun- 
dance ratios such [a/Fe] (where a-elements are O, Mg, Ne, Si, S, Ca and 
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Fig. 16. In this figure we show the smallest stellar mass which dies at any given 
[Fe/H] achieved by the ISM as a consequence of chemical evolution. Thus, it is clear 
that in the early phases of the halo only massive stars are dying and contributing 
to the chemical enrichment process. Clearly this graph depends upon the assumed 
stellar lifetimes and upon the age- [Fe/H] relation. It is worth noting that the Fe pro- 
duction from Type la SNe appears before the gas has reached [Fe/H] =-1.0, therefore 
during the halo and thick disk phase. This clearly depends upon the assumed Type 
la SN progenitors (in this case the single degenerate model). 



Ti) versus [Fe/H], a typical way of plotting the abundances measured in the 
stars. The time-delay refers to the delay with which Fe is ejected into the 
ISM by SNe la relative to the fast production of a-elements by core-collapse 
SNe. Tinsley (1979) first suggested that this time delay would have produced 
a typical signature in the [a/Fe] vs. [Fe/H] diagram. In the following years, 
Greggio & Renzini (1983b), by means of simple models (star formation burst 
or constant star formation) studied the effects of the delayed Fe production by 
Type la SNe on the [O/Fe] vs. [Fe/H] diagram. Matteucci & Greggio (1986) 
included for the first time the Type la SN rate formulated by Greggio & Ren- 
zini (1983a) in a detailed numerical model for the chemical evolution of the 
Milky Way. The effect of the delayed Fe production is to create an overabun- 
dance of O relative to Fe ([0/Fe]> 0) at low [Fe/H] values, and a continuous 
decline of the [O/Fe] ratio until the solar value ([O/Fe] = 0.0) is reached for 
[Fc/H]> —1.0 dex. This is what is observed and indicates that during the halo 
phase the [O/Fe] ratio is due only to the production of O and Fe by SNe II. 
However, since the bulk of Fe is produced by Type la SNe, when these latter 
start to be important then the [O/Fe] ratio begins to decline. This effect was 
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Fig. 17. The relation between [O/Fe] vs. [Fe/H] for Galactic stars in the solar 
vicinity. The models and the data are normalized to the solar meteoritic abundances 
of Anders & Grevesse (1989). The thick curve represents the predictions of the two- 
infall model where Type la SNe produce ~ 70% of Fe and Type II SNe the remaining 
~ 30%. The upper thin curve represents the case where all the Fe is assumed to be 
produced by Type la SNe, whereas the thin lower line refers to the case where all 
the Fe is assumed to be produced in Type II SNe. The data are from Melendez & 
Barbuy (2002). 

predicted by Matteucci & Greggio (1986) to occur also for other a-elements 
(e.g. Mg, Si). At the present time, a great amount of stellar abundances is 
available and the trend of the a-elements has been confirmed. Before showing 
some of the most recent data, it is worth showing better the time-delay model. 
In Figure 17 it is shown that a good fit of the [O/Fe] ratio as a function of 
[Fc/H] is obtained only if the a-elements are mainly produced by Type II 
SNe and the Fe by Type la SNe. If one assumes that only SNe la produce Fe 
as well as if one assumes that only Type II SNe produce Fe, the agreement 
with observations is lost. Therefore, the conclusion is that both Types of SNe 
should produce Fe in the proportions of 1/3 for Type II SNe and 2/3 for Type 
la SNe. The IMF also plays a role in this game and these proportions are ob- 
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taincd for "normal" Salpeter-like IMFs, which includes both Salpeter (1955) 
and Scalo (1986) or Kroupa et al. (1993) IMFs. 

As a further illustration of the time-delay model we show in Figures 18, 
19 and 20 the [X,/Fe] vs. [Fe/H] relations both observed and predicted for 
stars in the solar vicinity belonging to halo, thick- and thin-disk. The adopted 
yields for massive stars are those suggested by Frangois et al. (2004) in order 
to best fit these relations and the solar abundances (namely the abundances in 
the ISM 4.5 Gyr ago). These yields arc obtained by applying some corrections 
to the yields of WW95, as shown in Figure 21, where the ratios between the 
suggested and WW95 yields are reported. 

In Figure 22 we show the predictions of a chemical evolution model for 
the solar vicinity where the recent yields from massive stars of Nomoto et al. 
(2006) have been adopted. As one can sec, although some of the problems 
present in the previous yields have been alleviated, for other elements the 
disagreement still persists. 

The G-dwarf metallicity distribution and constraints on the thin 
disk formation 

The G-dwarf metallicity distribution is a quite important constraint for the 
chemical evolution of the solar vicinity. It is the fossil record of the star for- 
mation history of the thin disk. If one is able to reproduce such a distribution, 
then one can have an idea of the SFR and the IMF and, as a consequence, of 
the gas accretion history. Therefore, to fit the G-dwarf metallicity distribution 
means to obtain constraints on the mechanism of formation of the thin disk. 
Originally, there was the "G-dwarf problem" which means that the Simple 
Model of galactic chemical evolution could not reproduce the distribution of 
the G-dwarfs. It has been since long demonstrated that relaxing the closed-box 
assumption and allowing for the solar region to form gradually by accretion 
of gas can solve the problem (Tinsley, 1980; Pagel 1997). Also a variable IMF 
could solve the problem but it would create other problems (see Martinelli 
& Matteucci, 2000). Assuming that the disk forms from pre-enriched gas can 
also solve the problem but still the gas infall is necessary to have a realistic 
picture of the disk formation. The two-infall model can reproduce very well 
the G-dwarf distribution and also that of K-dwarfs (see Figures 23 and 24) , 
as long as a timescale for the formation of the disk in the solar vicinity of 
7-8 Gyr is assumed. This conclusion is shared by other authors (Alibes et al. 
2001; Prantzos & Boissier 1999) 

Carbon and Nitrogen evolution 

Carbon and nitrogen deserve a separate discussion from the other elements, 
in particular 14 N whose observational behaviour is difficult to reconcile with 
the theory. First of all, we should distinguish between primary and secondary 
elements: primary elements are those synthesized directly from H and He, 
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Fig. 18. Predicted and observed [a/Fe] vs. [Fe/H] in the solar neighbourhood. The 
models and the data are from Frangois et al. (2004). The models are normalized to 
the predicted solar abundances. The predicted abundance ratios at the time of the 
Sun formation (Solar value) are shown in each panel and indicate a good fit (all the 
values are close to zero). 



whereas secondary elements are those deriving from metals already present in 
the star at birth. In the framework of the Simple Model of galactic chemical 
evolution, the abundance of a secondary element evolves like the square of the 
abundance of the progenitor metal, whereas the evolution of the abundance 
of a primary element does not depend on the metallicity. 

In Figure 25 we show the predictions of the Simple Model for the ratio 
N/O, together with data for extragalactic HII regions and Damped Lyman-a 
systems (DLAs). 
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Fig. 19. The same as Fig. 18 for Ni, Zn, K and Sc. The models and the data are from 
Frangois et al. (2004). The models are normalized to the predicted solar abundances. 
The predicted abundance ratios at the time of the Sun formation are shown in each 
panel and indicate a good fit. 



It is worth noting that the solutions of the Simple Model for a primary 
and a secondary element are oversemplifications since the Simple Model does 
not take into account stellar lifetimes which are very important in 1A N pro- 
duction, which arises mainly from low and intermediate mass stars, both as 
a secondary and primary element (e.g. Renzini & Voli, 1981; van den Hoeck 
& Groenewegen 1997). Also 12 C originates mainly from low and intermediate 
mass stars. The contribution to 12 C from massive stars becomes very impor- 
tant only for metallicities oversolar, if the metallicity dependent mass loss is 
adopted (e.g. Maeder 1992). The interpretation of the diagram of Figure 25 is 
not so straightforward since extragalactic HII regions and DLAs are galaxies, 
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Fig. 20. The same as in Fig. 18 for Ti, Cr, Mn and Co. The models and the data 
are from Frangois et al. (2004). The models are normalized to the predicted solar 
abundances. The predicted abundance ratios at the time of the Sun formation are 
shown in each panel and indicate a good fit. 



and not necessarily that diagram is an evolutionary one, in the sense that 
O/H does not trace the time unlike [Fe/H] in the Galactic stars. Galaxies, 
in fact, may have started forming stars at different cosmic epochs and with 
different SF histories. However, if we interpret the diagram of Figure 25 as 
an evolutionary one, then the DLAs and the extragalactic HII regions of low 
metallicity should be younger and reflect the nucleosynthesis in massive stars 
and perhaps in intermediate mass stars. The observed plateau for N/O at low 
metallicity then would indicate a primary production of N in massive stars. 
Nitrogen, in fact, is also produced in massive stars: until a few years ago, the 
N production in massive stars was considered only a secondary process, until 
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Fig. 21. Ratios between the empirical yields derived by Frangois et al. (2004) and 
the yields of WW95 for massive stars. In the small panel at the bottom right we 
show the same ratios for SNe la and the comparison is with the yields of Iwamoto 
et al. (1999). 

Meynet & Maeder (2001, 2003, 2005) showed that stellar rotation in massive 
stars can produce primary N. A better test for the primary/secondary nature 
of N is represented by the Galactic stars, since they really represent an evo- 
lutionary sequence. In Figures 26 and 27 we show the most recent data on C 
and N compared with chemical evolution models including N from rotating 
massive stars. 

As one can see in Figures 26 and 27, the fit with data is good when 
primary N from massive stars is included. However, there are a few warnings, 
first of all the measurements of N abundance in stars of low metallicity are 
still uncertain and then the fact that the N measurement in the gas in DLAs 
at high redshift show that at low O abundances there are systems with a 
log(N/0) < —2.0, below the plateau shown by Galactic stars. A plateau in 
[N/Fe]is also observed in Galactic stars for [Fe/H]< —3.0 dex, as shown in 
Figure 27. In Figure 27 we show also the [C/Fe] values for Galactic stars but 
only for low metallicity stars: they indicate a roughly solar ratio like the stars 
with higher mctallicitics. Therefore, both [N/Fe] and [C/Fe] seem to show 
roughly constant solar values over the total [Fe/H] range. In the framework of 
the time-delay model, this means that C, N and Fe are all formed in the same 
stars (or with the same lifetimes) and that N is mainly a primary element. 
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Fig. 22. Predicted and observed [X/Fe] vs. [Fe/H] in the solar neighbourhood. The 
predictions are from Nomoto et al. (2006), where all the references to the data can 
be found, and they have been obtained by means of metal dependent yields. Figure 
from Nomoto et al. (2006). 



However, more data are necessary to assess this point and to reconcile the 
Galactic star data with high redshift DLAs. 

S- and r- process elements 

The s- and r- process elements are generally produced by neutron capture on 
Fe seed nuclei. The former are formed during the He-burning phase both in 
low and massive stars, whereas the latter occur in explosive events such as 
Type II SNe. Recently, Francois et al. (2006) have measured the abundances of 
several very heavy elements (e.g. Ba and Eu) in extremely metal poor stars of 
the Milky Way. Previous work on the subject had shown a large spread in the 
abundance ratios of these elements to iron, especially at low metallicities. This 
spread is confirmed by this more recent study although is less than before, 
and is at variance with the lack of spread observed in the other elements 
shown before (e.g. a-elements). Apart from this problem, not yet solved, these 
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Fig. 23. The G-dwarf metallicity distribution. The model prediction is from Chi- 
appini et al. (1997) and assumes a timescale for the formation of the local disk of 8 
Gyr. The data are represented by the histograms. 



diagrams can be very useful to place constraints on the nucleosynthetic origin 
of these elements. In particular, Cescutti et al. (2006) by adopting the two- 
infall model predicted the evolution of [Ba/Fe] and [Eu/Fe] versus [Fe/H], as 
shown in Figures 28 and 29. They can well fit the average trend but not the 
spread at very low mctallicities since the model assumes instantaneous mixing. 
In order to fit the Ba evolution, they assumed that Ba is mainly produced as 
s- process element in low mass stars (1-3M Q ) but that a fraction of Ba is also 
produced as an r-process element in stars with masses 12-30Mq. Europium is 
assumed to be only an r-process element produced in the range 12-30M Q . 

In order to explain why the s- and r- process elements show a large and 
probably real spread at very low metallicities, whereas elements such as the 
a-elements show only a little spread, one could think of a moderately inho- 
mogeneous model coupled with differences in the nucleosynthesis between s- 
and r- process elements on one side and cv-elements on the other side (see 
Cescutti 2008). Highly inhomogeneous models for the halo evolution, in fact, 
predict a too large spread for the a-elements at low metallicity (e.g. Argast 
et al. 2000). It is worth noting the typical secondary behaviour of Ba, whose 
main production is by means of the s-process, which needs Fe seed nuclei al- 
ready present in the star, and neutrons which are accreted on these nuclei. 
The production of neutrons is also dependent on the original metal content, 
therefore it would be even more precise to speak of Ba as a tertiary element. 
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Fig. 24. The figure is from Kotoneva et al. (2002) and shows the comparison be- 
tween a sample of K-dwarfs and model predictions in the solar neighbourhood. The 
dotted curve refers to the two-infall model with a timescale r = 2 Gyr, whereas the 
continuous line refers to r = 8 Gyr, as in Fig. 23. 



3.5 The Galactic disk 

A good model of chemical evolution for the Milky Way should reproduce also 
the features of the Galactic disk. In particular: abundance gradients, gas and 
SFR distribution with the galactocentric distance. 

Abundance gradients 

The chemical abundances measured along the disk of the Galaxy suggest that 
the metal content decreases from the innermost to the outermost regions, in 
other words there is a negative gradient in metals. Abundance gradients can 
be derived from HII regions, planetary nebulae (PNe), open clusters and stars 
(0,B stars and Cepheids). There are two types of abundance determinations 
in HII regions: one is based on recombination lines which should have a weak 
dependence on the temperature of the nebula (He, C, N, O), the other is 
based on collisionally excited lines where a strong dependence is intrinsic to 
the method (C, N, O, Ne, Si, S, CI, Ar, Fe and Ni). This second method has 
predominated until now. A direct determination of the abundance gradients 
from HII regions in the Galaxy from optical lines is difficult because of ex- 
tinction, so usually the abundances for distances larger than 3 Kpc from the 
Sun arc obtained from radio and infrared emission lines. 
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Fig. 25. The plot of log (N/O) vs. log(0/H) + 12: small dots represent extragalactic 
HII regions, red triangles are Damped-Lyman a systems (DLA), which are high 
redshift objects (Pettini et al. 2002). The three green triangles are the most recent 
determinations for DLAs from Pettini et al. (2008). Dashed lines mark the solution 
of the Simple Model for a primary and a secondary element. Figure from Pettini et 
al. (2008). 
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Fig. 26. Upper panel: solar vicinity diagram log(N/0) vs. log(0/H) + 12. The data 
points are from Israelian et al. (2004) (large squares) and Spite et al. (2005) (aster- 
isks). Models: the dashed line represents a model with substantial primary N pro- 
duction from massive stars. This was obtained by means of stellar models (Meynet et 
al. 2006; Hirschi 2007) with faster rotation relative to the work of Meynet & Maeder 
(2002) for Z = 10~ 8 . Lower panel: solar vicinity diagram log(C/0) vs. log(0/H) + 
12. The data are from Spite et al. (2005) (asterisks), Israelian et al. (2004) (squares) 
and Nissen (2004) (filled pentagons). Solar abundances (Asplund et al. 2005, and 
references therein) are also shown. Figure from Chiappini et al. (2006). 
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Fig. 27. Observed and predicted [C/Fe] vs. [Fe/H] (upper panel) and [N/Fe] vs. 
[Fe/H] (bottom panel) in the solar neighbourhood. The data points are from Cayrel 
et al. (2004), Spite et al. (2005) (asterisks) and Israelian et al. (2004) (squares). The 
dot-dashed line represents a model with yields from Chiefh & Limongi (2002, 2004) 
for a metallicity Z = 10~ 6 connected to the Pop III stars (only massive stars for 
that metallicity). The dashed line and the dotted lines represent heuristic models 
where the yields of C and N have been assumed "ad hoc" . In particular, the fraction 
of primary N from massive stars is obtained by the fit to the data at low metallicity. 
Figure from Chiappini, Matteucci & Ballero (2005). 
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Fig. 28. The evolution of Barium in the solar vicinity as predicted by the two-infall 
model (Cescutti et al. 2006). Data are from Francois et al. (2006). 

Abundance gradients can also be derived from optical emission lines in 
PNe. However, the abundances of He, C and N in PNe are giving only infor- 
mation on the internal nucleosynthesis of the star. So, to derive gradients one 
should look at the abundances of O, S and Ne, unaffected by stellar processes. 
Abundance gradients are derived also from measuring the Fe abundance in 
open clusters (e.g. Carraro et al. 2004; Yong et al. 2005) or from abundances 
in Cepheids (e.g. Andrievsky et al. 2002 a,b,c, 2004; Luck et al. 2003; Yong 
et al. 2006) or from abundances in O, B stars (e.g. Daflon & Cunha, 2004). 

In Figure 30 we show theoretical predictions of abundance gradients along 
the disk of the Milky Way compared with data from HII regions, B stars and 
PNe. The adopted model is from Chiappini et al. (2001) and is based on an 
inside-out formation of the thin disk. The assumed model does not allow for 
exchange of gas between different regions of the disk. The disk is, in fact, 
divided in several concentric shells 2 Kpc wide with no interaction between 
them. 

As already mentioned, most of the current models agree on the inside-out 
scenario for the disk formation, however not all models agree on the evolution 




Fig. 29. The evolution of Europium in the solar vicinity (Cescutti et al. 2006). Data 
are from Frangois et al. (2006). 



of the gradients with time. In fact, some models, although assuming an inside- 
out formation of the disk, predict a gradient flattening with time (Boissier & 
Prantzos 1999; Alibes et al. 2001), whereas others such as that of Chiappini 
et al. (2001) predict a steepening, as shown in Figure 30. The reason for the 
steepening is that in the model of Chiappini et al. there is included a threshold 
density for SF, which induces the SF to stop when the density decreases 
below the threshold. This effect is particularly strong in the external regions 
of the Galactic disk, thus contributing to a slower evolution and therefore to a 
steepening of the gradients with time. In Figure 31 we show models and some 
more recent data including Cepheids. 

In the Chiappini et al. model, the fit to the gradients is obtained by 
means of the inside-out formation of the Galactic disk. Numerical simula- 
tions of abundance gradients show that no gradient arises if one assumes the 
same timescale of disk formation at any galactocentric distance. The differ- 
ent timescale of accretion influences the SFR, thus creating a gradient in the 
SFR and therefore in the resulting metal content. However, it should be said 
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Fig. 30. Spatial and temporal behaviour of abundance gradients along the Galactic 
disk as predicted by the best model of Chiappini et al. (2001). The upper lines in 
each panel represent the present time gradient, whereas the lower ones represent the 
gradient a few Gyr ago. It is clear that the gradients tend to stepeen in time, a still 
controversial result. The data are from HII regions, B stars and PNe (see Chiappini 
et al. 2001). 

that the effect of the threshold is also important and tends to steepen the 
gradients. 

In Figure 32 we show the results of Boissier & Prantzos (1999) for abun- 
dance gradients and also for the gas and SFR distribution along the disk. Note 
the gradients are flattening with time. 

3.6 The Galactic bulge 
Bulge formation 

The bulges of spiral galaxies are generally distinguished in true bulges, hosted 
by SO-Sb galaxies and "pseudobulges" hosted in later type galaxies (see Ren- 
zini 2006 for references). Generally, the properties (luminosity, colors, line 
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Fig. 31. Gradients of the a-elements along the disk. The predicted gradients for 
O, Mg, Si, S and Ca are compared with different sets of data. The small open circles 
are the data of the Cepheids by Andrievsky et al. (2002a,b,c; 2004) and Luck et al. 
(2003). The solid triangles are the data by Daflon & Cunha (2004) (OB stars), the 
open squares are the data by Carney et al. (2005) (red giants), the solid hexagons 
are the data by Yong et al. (2006) (Cepheids), the open triangles are the data by 
Yong et al. (2005) (open clusters) and the solid squares are the data by Carraro et 
al. (2004) (open clusters) . The most distant value for Carraro et al. (2004) and Yong 
et al. (2005) refers to the same object: the open cluster Berkeley 29. The thin solid 
line represents the model predictions at the present time normalized to the mean 
value of the Cepheids at 8 Kpc; the dashed line represents the predictions of the 
model at the epoch of the formation of the solar system normalized to the observed 
solar abundances by Asplund et al. (2005). This prediction should be compared with 
the data for red giant stars and open clusters (Carraro et al. 2004; Carney et al. 
2005; Yong et al. 2005). The models and the Figure are from Cescutti et al. (2007). 
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Fig. 32. Comparison between model predictions and observations for the disk of 
the Milky Way. The figure is from Boissier & Prantzos (1999). Top left panel: gas 
distribution along the disk. Top right panel: the O gradient at the present time (curve 
with label 13.5) and at two other different cosmic epochs (5 Gyr and 1 Gyr from the 
beginning) . Second left panel: the surface mass density of living stars. Second right 
panel: the Fe gradient. Third left panel: the gradient of the SFR normalized to the 
value at the solar ring. Third right panel: the predicted distribution of the current 
surface mass densities of stellar remnants (WDs) , black holes (BH) and neutron stars 
(NS). Fourth left panel: the predicted infall rate along the disk at three different 
cosmic epochs. Fourth right panel: the predicted distributions of surface densities 
by number of the stellar remnants. 
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strenghts) of true bulges are very similar to those of elliptical galaxies. In the 
following, we will refer only to true bulges and in particular to the bulge of 
the Milky Way. The bulge of the Milky Way is, in fact, the best studied bulge 
and several scenarios for its formation have been put forward in past years. 
As summarized by Wyse & Gilmore (1992) the proposed scenarios are: 

• the bulge formed by accretion of extant stellar systems which eventually 
settle in the center of the Galaxy. 

• The bulge was formed by accumulation of gas at the center of the Galaxy 
and subsequent evolution with either fast or slow star formation. 

• The bulge was formed by accumulation of metal enriched gas from the 
halo, thick disk or thin disk in the Galaxy center. 

In the context of chemical evolution, the Galactic bulge was first modeled 
by Matteucci & Brocato (1990) who predicted that the [a/Fe] ratio for some 
elements (O, Si and Mg) should be supersolar over almost the whole metallic- 
ity range, in analogy with the halo stars, as a consequence of assuming a fast 
bulge evolution which involved rapid gas enrichment in Fe mainly by Type II 
SNc. At that time, no data were available for chemical abundances; the pre- 
dictions of Matteucci & Brocato (1990) were confirmed for a few a-elements 
(Mg, Ti) by the observations of McWilliam & Rich (1994, hereafter MR94), 
whereas for other a-elements (e.g. Ca, Si) the observed trend was different. 
Other discrepancies regarding the Mg overabundance came from Sadler et al. 
(1996). In order to better assess these points, Matteucci et al. (1999) studied 
a larger set of abundance ratios, by means of a detailed chemical evolution 
model whose parameters were calibrated so that the metallicity distribution 
observed by MR94 could be fitted. They concluded that an evolution much 
faster than that in the solar neighbourhood and even faster than that of the 
halo (see also Renzini, 1993) is necessary for the MR94 metallicity distribu- 
tion to be reproduced, and that an IMF index flatter (x = 1.1 — 1.35) than 
that of the solar neighbourhood is needed as well. They also made predic- 
tions about the evolution of several abundance ratios which were meant to be 
confirmed or disproved by subsequent observations, namely that a-elemcnts 
should in general be overabundant with respect to Fe, but some (e.g. Si, Ca) 
less than others (e.g. O, Mg), and that the [ 12 C/Fe] ratio should be solar at 
all metallicities. 

Samland et al. (1997) developed a self-consistent chemo-dynamical model 
for the evolution of the Milky Way components starting from a rotating pro- 
togalactic gas cloud in virial equilibrium, which collapses owing to dissipative 
cloud-cloud collisions. They found that self-regulation due to a bursting star 
formation and subsequent injection of energy from Type II supernovae led 
to the development of "contrary flows", i.e. alternate collapse and outflow 
episodes in the bulge. This caused a prolonged star formation episode lasting 
over ~4x 10 9 yr. They included stellar nucleosynthesis of O, N and Fe, but 
claimed that gas outflows prevent any clear correlation between local star for- 
mation rate and chemical enrichment. With their model, they could reproduce 



46 Francesca Matteucci 



the oxygen gradient of HII regions in the equatorial plane of the Galactic disk 
and the metallicity distribution of K giants in the bulge (Rich, 1988), field 
stars in the halo and G dwarfs in the disk, but they did not make predic- 
tions about abundance ratios in the bulge. In general, hierarchical clustering 
models of galaxy formation do not support the conclusion of a fast formation 
and evolution of the bulge. In Kauffmann (1996) the bulges form through 
violent relaxation and destruction of disks in major mergers. The stars of the 
destroyed disk build the bulge, and subsequently the bulge has to be rebuilt. 
This implies that late type spirals should have older bulges than early type 
ones, since the build-up of a large disk needs a long time during which the 
galaxy has to evolve undisturbed. This is not confirmed by observations, since 
the high metallicity and the the narrow age distribution observed in bulges of 
local spirals are not compatible with their merger origin (see Wyse, 1999). 

Molla et al. (2000) proposed a multiphase model in the context of the 
dissipative collapse scenario of the Eggen et al. (1962) picture. They supposed 
that the bulge formation occurred in two main infall episodes, the first from 
the halo to the bulge, on a timescale th = 0.7 Gyr, and the second from 
the bulge to a so-called core population in the very nuclear region of the 
Galaxy, on a timescale tb 3> th- The three zones (halo, bulge, core) interact 
via supernova winds and gas infall. They concluded that there is no need for 
accretion of external material to reproduce the main properties of bulges and 
that the analogy to ellipticals is not justified. Because of their rather long 
timescale for the bulge formation, these authors did not predict a noticeable 
difference in the trend of the [a/Fe] ratios but rather suggested that they 
behave more akin to that in the solar neighbourhood (contrary to several 
indications from abundance data, e.g. MR94). 

Immeli et al. (2004) investigated the role of cloud dissipation in the forma- 
tion and dynamical evolution of star forming gas rich disks by means of a 3D 
chemodynamical model. They found that the galaxy evolution proceeds very 
differently depending on whether the gas disk or the stellar disk first become 
unstable. This in turn depends on how efficiently the cold cloud medium can 
dissipate energy. If the gas cools efficiently, a starburst takes place which gives 
rise to enhanced [a/Fe] ratios, thus in agreement with a fast bulge formation. 

A more recent model was proposed by Costa et al. (2005), in which the 
best fit to observations is achieved by means of a double infall model. An 
initial fast (0.1 Gyr) collapse of primordial gas is followed by a supernova- 
driven mass loss and then by a second slower (2 Gyr) infall episode, enriched 
by the material ejected by the bulge during the first collapse. Costa et al. 
(2005) claimed that the mass loss is necessary to reproduce the abundance 
distribution observed in PNe, and because the predicted abundances would 
otherwise be higher than observed. However, it should be noted again that 
the abundances derived from PNe can be affected by internal stellar processes 
and therefore are meaningless for studying galactic chemical evolution. With 
their model, they are able to reproduce the trend of [O/Fe] abundance ratio 
observed by Pompeia et al. (2003) and the data of nitrogen versus oxygen 
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abundance observed by Escudero & Costa (2001) and Escudero et al. (2004). 
It must be noted however that Pompeia et al. (2003) obtained abundances for 
"bulge-like" dwarf stars. This "bulge-like" population consists of old (~ 10—11 
Gyr) metal-rich nearby stars whose kinematics and metallicity suggest an 
inner disk or bulge origin and a mechanism of radial migration, perhaps caused 
by the action of a Galactic bar, but the birthplace of these stars is undoubtedly 
not certain. 

Interpretation of bulge data and other galaxies 

In summary, MR94 first measured the metallicity distribution and the [a/Fe] 
ratios in the Galactic bulge and confirmed partly the predictions of Mat- 
teucci & Brocato (1990) that all of the a-elements should be enhanced rela- 
tive to Fe for a large range of [Fc/H]. In fact, MR94 found that not all the 
a-elements were enhanced, in particular oxygen. Very recently, medium- and 
high-resolution spectroscopy of bulge stars was performed (Rich & Mc William, 
2000; Fulbright et al., 2006, 2007; Zoccali et al., 2006; Lecureur et al. 2007), 
and it seems to indicate that also O is enhanced, thus supporting the sug- 
gestion of a fast formation of the bulge. The metallicity distribution of stars 
in the bulge and the [a/Fe] ratios greatly help in selecting the most proba- 
ble scenario for the bulge formation. In Figure 33 we present the predictions 
by Matteucci (2003) of the [a/Fe] ratios as functions of [Fe/H] in galaxies 
of different morphological type. In particular, for the Galactic bulge or an 
elliptical galaxy of the same mass, for the solar vicinity region and for an 
irregular magellanic galaxy (LMC and SMC). The underlying assumption is 
that different objects undergo different histories of star formation, being very 
fast in the spheroids (bulges and ellipticals), moderate in spiral disks and slow 
and perhaps gasping in irregular gas rich galaxies. The effect of different star 
formation histories is evident in Figure 33 where the predicted [a/Fe] ratios in 
the bulge and ellipticals remain high and almost constant for a large interval 
of [Fe/H]. This is due to the fact that, since star formation is very intense, the 
bulge reaches very soon a solar metallicity thanks only to the SNe II; then, 
when SNe la start exploding and ejecting Fe into the ISM, the change in the 
slope occurs at larger [Fe/H] than in the solar vicinity. In the extreme case of 
irregular galaxies the situation is opposite: here the star formation is slow and 
when the SNe la start exploding the gas is still very metal poor. This scheme 
is quite useful since it can be used to identify galaxies only by looking at their 
abundance ratios. A model for the bulge behaving as shown in Figure 33 is 
able to reproduce also the observed metallicity distribution of bulge stars (see 
Matteucci & Brocato 1990; Matteucci et al. 1999). The scenario suggested in 
these papers favors the formation of the bulge by means of a short and strong 
starburst, in agreement with Elmegreen (1999) and Ferreras et al. (2003). A 
similar model, although updated with the inclusion of the development of a 
galactic wind and more recent stellar yields, has been presented by Ballero 
et al. (2007): it shows how a model with intense star formation (star for- 
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mation efficiency ~ 20Gyr~ 1 ) and rapid assembly of gas (0.1 Gyr) can best 
reproduce the most recent accurate data on abundance ratios and metallicity 
distribution. 
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Fig. 33. The predicted [a/Fe] vs. [Fe/H] relations for the Galactic bulge (upper 
curve), the solar vicinity (median curve) and irregular galaxies (low curve). Data for 
the bulge are reported for comparison. Data for the LMC and DLA systems are also 
shown for comparison, indicating that DLAs are probably irregular galaxies. Figure 
and references are from Matteucci (2003). 



In Figures 34 and 35 we show a detailed comparison between model pre- 
dictions for the Galactic bulge and data on O and Mg. 

As one can see, the plateau in the [a/Fe] is longer than in the solar neigh- 
bourhood, since in the bulge the slope of the [a/Fe] ratio starts changing 
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Fig. 34. The predicted [O/Fe] vs. [Fe/H] relation for the Bulge (curve), compared 
with the most recent data. The chemical evolution model is that of Ballero et al. 
(2007), where references to the data can be found. 



drastically only for [Fe/H]> 0.0 dex. The long plateau is well explained by 
the model of Ballero et al. (2007) assuming a very fast formation of the bulge. 
It is worth noting that the [O/Fe] ratio has a steeper slope than [Mg/Fe] and 
this could be due to differences in the nucleosynthesis of these elements (e.g. 
McWilliam et al. 2007). 

The IMF assumed for the bulge is usually flatter than the IMF of the solar 
neighbourhood and this is generally dictated by the fit of the bulge metallicity 
distribution which peaks at a higher [Fe/H] than the G-dwarf metallicity 
distribution in the solar vicinity. Numerical calculations have indicated that 
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Fig. 35. The predicted [Mg/Fe] vs. [Fe/H] relation for the bulge (curve), compared 
with the most recent data. The chemical evolution model is that of Ballero et al. 
(2007), where references to the data can be found. 



the main parameter influencing the peak of the distribution is the IMF, as 
clearly shown in Figures 36 and 37. 

In summary, the comparison between the models (Ballero et al. 2007) 
on one side, and the metallicity distribution and the [a/Fe] ratios on the 
other, strongly indicates that the Galactic bulge is very old and must have 
formed very quickly during a strong starburst (with a SF efficiency much 
higher than in the disk). The metallicity distribution in particular, seems to 
suggest an IMF flatter than in the disk with an exponent for massive stars 
in the range x = 1.35 — 0.95. However, to assess more precisely this point we 
need more data: in particular, a flatter IMF predicts that the overabundances 
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Fig. 36. The predicted and observed metallicity distribution in the Galactic bulge. 
The data are from Zoccali et al. (2003) (dashed histogram) and Fulbright et al. 
(2006) (continuous histogram). In particular, the model with the peak at the lower 
metallicity is computed with an IMF which is similar to that of the solar vicinity and 
indicated by IGIMF, whereas the distribution which best fits the data is computed 
with a flat IMF (x=0.95 for M > 1M©). 

of a-elements relative to Fe and to the Sun should be higher in the bulge 
than in the disk. This is not entirely clear from the available data, although 
Zoccali et al. (2006) conclude that the [O/Fe] ratios in bulge stars are higher 
than in thick and thin disk stars (see Figure 38). The timescale for the bulge 
formation by accretion of gas lost from the halo is 0.1 Gyr and certainly no 
longer than 0.5 Gyr. 

4 What we have learned about the Milky Way 

From the discussions of the previous sections we can extract some important 
conclusions on the formation and evolution of the Milky Way, derived from 
chemical abundances. In particular: 

• The inner halo formed on a timescale of 1-2 Gyr at maximum, the outer 
halo formed on longer timcscales perhaps from accretion of satellites or 
gas. 

• The disk at the solar ring formed on a timescale not shorter than 7 Gyr. 
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Fig. 37. The predicted and observed metallicity distribution in the Galactic bulge. 
The data are from Zoccali et al (2003) (dashed histogram) and Fulbright et al. (2006) 
(continuous histogram). The lines are the predictions of models with the same IMF 
but different SF efficiencies, as indicated in the Figure. 



• The whole disk formed inside out with timescales of the order of 2 Gyr or 
less in the inner regions and 10 Gyr or more in the outermost regions. 

• The abundance gradients arise naturally from the assumption of the inside- 
out formation of the disk. A threshold density for the star formation helps 
in steepening the gradients in the outer disk regions. 

• The bulge is very old and formed very quickly on a timescale smaller than 
even the inner halo and not larger than 0.5 Gyr. 

• The IMF seems to be different in the bulge and the disk, being flatter in 
the bulge, although more abundance data are necessary before drawing 
firm conclusions. 



5 The time-delay model and the Hubble sequence 

In this section we will discuss how different star formation histories affect the 
evolution of galaxies of different morphological type and in particular how the 
abundance patterns are expected to change with the star formation. 
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Fig. 38. Upper panel: the [O/Fe] vs. [Fe/H] for bulge stars as measured by Zoccali 
et al. (2006) (circles), along with previous determinations in other bulge stars from 
optical (open triangles) and near-IR spectra (crosses). Open symbols refer to spec- 
tra with lower S/N or with the O line partially blended with telluric absorption. 
Lower panel: [O/Fe] trend in the bulge and in thick and thin disk stars (crosses and 
triangles). The solid line shows a linear fit to the thick data points with [Fe/H]> 
-0.5 dex and is meant to emphasize that all bulge stars with— 0.4 < [Fe/H] < +0.1 
are more O-enhanced than the thick disk stars. If confirmed, this trend may indicate 
a flatter IMF for the bulge than for the disk. This plot also suggests that there is a 
systematic difference between bulge and disk stars, thus excluding that bulge stars 
were once disk stars migrated in the bulge. Figure and references from Zoccali et al. 
(2006). 
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5.1 Star formation and Hubble sequence 

Sandage (1986), on the basis of a work of Gallagher, Hunter & Tutukov (1984) 
who measured SFR in galaxies, suggested a possible interpretation of the 
Hubble sequence in terms of different star formation histories. In this picture, 
ellipticals and bulges must have suffered an intense and strongly declining 
SFR, whereas late type galaxies must have undergone through a less intense, 
almost constant (spirals) and even increasing with time (irregulars) SFR. In 
Figure 39 is illustrated such a behaviour of the SFR for E, SO and Sa galaxies: 
there are two important timescales t c , the gas collapse time, and t s , the star 
formation time, namely the timescale on which the gas in a galaxy is consumed 
by means of star formation (the inverse of the SF efficiency). The interplay 
between these two quantities can be crucial for the formation of the different 
galactic morphological types. In fact, if t s << t c , most of the stars form before 
the collapse is over and the gas does not have time to dissipate energy and 
settle into a disk. In this case, the resulting galaxy will be a spheroid, whereas 
if t s > t c the gas has time to dissipate energy and form a spiral galaxy. 
This picture is certainly too simplistic to fully describe the reality of galaxy 
formation but it seems to work well when we are going to interpret galaxy 
formation by studying the stellar populations in galaxies. 

Later on, Kennicutt (1998a, b) measured again the SFR in star forming 
galaxies and suggested similar behaviours for the different galactic morpho- 
logical types, as shown in Figure 40. In this Figure, the behaviour of the SFR 
in spirals is obtained by fitting the mean value of the parameter 6, measured 
by Kennicutt et al. (1994). The parameter 6 is the ratio between the present 
time SFR and the average SFR in the past, as defined in eq.(16). Early type 
spiral galaxies are characterized by rapidly declining SFRs, with b ~ 0.01 — 0.1, 
whereas late type spirals have formed stars since a long time at an almost con- 
stant rate with 6=1. Finally, ellipticals and SO have long ago ceased forming 
stars and have 6 = 0. 

In Figure 41 we show the SF histories which give rise to the [a/Fe] vs. 
[Fe/H] relations of Fig. 33. 

In Figure 42 we show the predicted Type la SN rates according to the 
SFRs of Figure 41. The assumed progenitor model for Type la SNe is the 
single degenerate with the delay time distribution as in Matteucci & Recchi 
(2001). 

The typical timescale for SN la enrichment 

The predicted Type la SN rates for galaxies with different morphologies show a 
difference in the maximum SN la rate which is reached quite early in ellipticals 
and it occurs later and later moving to late types. Matteucci & Recchi (2001) 
suggested to assume the time for the occurrence of the maximum SNIa rate as 
the typical timescale for the chemical enrichment from these SNe. It depends 
on the star formation history of each galaxy, on the IMF and on the stellar 
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Fig. 39. Schematic representation of the change of the SFR with time for galaxies of 
types E, SO, and Sa. The dashed vertical line at the collapse time t c separates regions 
of low energy dissipation (to the left) from those of high energy dissipation (to the 
right). Bulges form in the left region, disks in the right. The integral under the 
curves gives the total number of stars formed (per unit galaxy mass). The integral 
under the SO curve is shaded. Figure from Sandage (1986). 



lifetimes. As we have already shown, the IMF together with stellar lifetimes 
represent the distribution of the time-delays (DTD) with which SNe explode; 
therefore, when a DTD is specified, the Type la SN rate depends on the SF 
history only. In summary, in ellipticals and bulges this timescale is 0.3-0.5 Gyr 
since the beginning of star formation, in the solar vicinity there is a first peak 
at 1 Gyr then it decreases slightly and increases again till 3 Gyr (due to the 
two-infall episodes). In irregular galaxies the maximum is reached at ~ 4 Gyr 
and then the rate remains constant. 
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Fig. 40. Schematic illustration of the stellar birthrate for different Hubble types. 
The left panel shows the evolution of the relative SFR with time following Sandage 
(1986). The curves for spiral galaxies are exponentially declining SFRs which fit the 
mean values of the birthrate parameter b as measured by Kennicutt et al. (1994). 
The curve for elliptical galaxies and bulges is an arbitrary dependence for an e- 
folding time of 0.5 Gyr, for comparative purposes only. The right panel shows the 
corresponding evolution in SFR with redshift, for an assumed cosmological density 
parameter J? = 0.3 and formation redshift Zf = 5. Figure from Kennicutt (1998b). 



6 Dwarf Spheroidals of the Local Group 

A different pattern for the [a/Fe] vs. [Fe/H] relation compared to the solar 
vicinity is observed in dwarf spheroidal galaxies (dSphs) of the Local Group, 
as shown in Figure 43, and this can be easily interpreted in the framework of 
the time-delay model coupled with different star formation histories. 

Before interpreting the [a/Fe] diagram, we recall the current ideas about 
the formation of the dSphs. 




Fig. 41. The predicted histories of SF in galaxies of different morphological type, 
with decreasing efficiency of SFR from ellipticals to irregulars. Figure from Calura 
(2004). 



6.1 How do dSphs form? 

Cold dark matter (CDM) models for galaxy formation predict that the dSphs, 
systems with luminous masses of the order of 10 7 Mq, are the first objects to 
form stars and that all stars in these systems should form on a timescale < 1 
Gyr, since the heating and gas loss, due to reionization, must have halted 
the SF soon. However, observationally all dSph satellites of the Milky Way 
contain old stars indistinguishable from those of Galactic globular clusters 
and they seem to have experienced SF for long periods (> 2 Gyr, Grebel & 
Gallagher 2004). The histories of SF for these galaxies are generally derived 
from the color-magnitude diagram (e. g. Mateo 1998). By looking at the [a/Fe] 
vs. [Fe/H] relations for dSphs, as shown in Figure 43, one can immediately 
suggest, on the basis of the time-delay model, that their evolution should 
have been characterized by a slow and protracted SF, at variance with the 
suggestion of a fast episode of SF truncated by the heating due to reionization. 
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Fig. 42. Predicted Type la SN rates (expressed in Gyr" 1 obtained according to 
the SF histories of Figure 41. As one can see the elliptical (dotted line) reaches 
a maximum at 0.3 Gyr, whereas the spiral (continuous line) at ~ 2 Gyr and the 
irregular (dashed line) has a rate increasing to its maximum at around 5 Gyr and 
then is roughly constant for the rest of the galactic lifetime. The minimum time delay 
for the Type la SNe to appear is 30 Myr, hardly visible in this plot. 
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Fig. 43. Observed [a/Fe] vs. [Fe/H] in the Milky Way (small points) and in dSphs 
(points with error bars). Figure from Shetrone et al. (2001). 



Dark matter in dSphs 

The dSph satellites of the Milky Way are considered the smallest dark matter 
dominated systems in the universe. In the past years there have been a few 
attempts at deriving the amount of dark matter in dSphs, in particular by 
measuring the mass to light ratios versus magnitude for these galaxies (e.g. 
Mateo 1998; Gilmore et al. 2007). Gilmore et al. (2007) suggested that the 
dSphs have a shallow central dark matter distribution and no galaxy is found 
with a dark mass halo less massive than 5 • 10 7 M Q , as shown in Figure 44. 



Observations of dSphs 

In recent years there has been a fast development in the field of chemical 
evolution of dSphs of the Local Group due to the increasing amount of data 
on chemical abundances derived from high resolution spectra (e.g. Smccker- 
Hane & McWilliam, 1999; Bonifacio el al. 2000, 2004; Shetrone et al. 2001, 
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Fig. 44. Dark matter in dSphs: mass to light ratios versus absolute V magnitude 
for some Local Group dSphs. The solid curve shows the relation expected if all the 
dSphs contain about 4 • 10 7 Mq of dark matter interior to their stellar distributions. 
Figure from Gilmore et al. (2007). 



2003; Tolstoy et al.2003; Bonifacio et al. 2004; Venn et al.2004; Sadakane et 
al. 2004; Fulbright et al. 2004; McWilliam & Smecker-Hane 2005a,b; Monaco 
et al. 2005; Geisler et al. 2005). The abundances of a-elements (O, Mg, Ca, 
Si,) plus the abundances of s- and r- process elements (Ba, Y, Sr, La and Eu) 
were measured with unprecedented accuracy. Besides these high-resolution 
studies we recall also the measure of the metallicities of many red giant stars 
in several dSphs by Tolstoy et al. (2004), Koch et al. (2006), Helmi et al. 
(2006) and Battaglia et al. (2006) obtained from the low-resolution Ca triplet. 
An interesting result of Helmi et al. (2006) is that they did not find stars 
with [Fe/H]< -3.0 dex and that the metallicity distribution of the stars in 
dSphs is different from that of halo stars in the Milky Way. Other important 
information, as already mentioned, comes from the photometry of dSphs of the 
Local Group and in particular from the color-magnitude diagrams. From these 
diagrams one can infer the history of SF of these galaxies. We recall the studies 
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of Hernandez et al. (2000), Dolphin (2002), Bellazzini ct al. (2002), Rizzi et 
al. (2003), Monelli et al. (2003), Dolphin ct al. (2005). The color-magnitude 
diagrams seem to indicate that the majority of dSphs had one rather long 
episode of SF with the exception of Carina for which four episodes of SF have 
been suggested (Rizzi et al. 2003). 

Chemical evolution of dSphs 

Several papers have appeared in the last few years concerning the chemical 
evolution of dSphs. For example Carigi et al. (2002) computed models for the 
chemical evolution of four dSphs by adopting the SF histories derived, from 
color-magnitude diagrams, by Hernandez ct al. (2000). In their model they 
assumed gas infall and computed the gas thermal energy heated by SNe in 
order to study galactic winds. In fact, the dSphs must have lost their gas in 
one way or another (galactic winds and/or ram pressure stripping) since they 
appear completely without gas. They assumed that the wind is sudden and 
devoids the galaxy of gas instantaneously. The adopted IMF is the Kroupa et 
al. (1993) IMF as in the solar vicinity. They predicted a too high metallicity 
for Ursa Minor and did not match the correct slope for the observed [a/Fe] 
ratios, unless the history of SF in this galaxy was assumed to be different than 
suggested by its color- magnitude diagram, as shown in Figure 45. 

Then, Ikuta & Arimoto (2002) proposed a closed box model (no infall nor 
outflow) for dSphs. In this Simple Model they had to assume some external 
cause to stop star formation, such as ram pressure stripping. They tested 
different IMFs and suggested that these galaxies had suffered very low SFRs 
(1-5% of that in the solar neighbourhood) and that the SF had a long duration 
(> 3.9-6.5Gyr). In Figure 46 are shown their predictions for [Mg/Fe] in dSphs. 
Also here, the predicted slope of the [Mg/Fe] ratio is flatter than observed. 

More recently, Fenner et al. (2006) suggested a model with a galactic wind 
for Sculptor: they indicated an efficiency of SF of O.ObGyr^ 1 . They concluded, 
from the study of the [Ba/Y] ratio, that chemical evolution in dSphs is in- 
consistent with the SF being truncated after reionization (at redshift z=8). In 
fact, the high value of this ratio measured in stars indicates strong s-process 
production from low mass stars which have very long lifetimes. 

The results of Lanfranchi &: Matteucci 

Lanfranchi & Matteucci (2003; 2004 hereafter LM04) developed models for 
dSphs of the Local Group. First they tested a "standard model" devised for 
describing an average dSph galaxy. This model was based on the following 
assumptions: 

• one long star formation episode of duration ~ 8 Gyr, 

• a small star formation efficiency, namely the star formation rate per unit 
mass is 1-10% of that in the solar vicinity, 
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Fig. 45. Observed and predicted [O/Fe] vs. [Fe/H] relation for the galaxy Ursa 
Minor. Models and Figures by Carigi et al. (2002). The upper Figure presents models 
assuming only one burst of SF at high redshift, as suggested by the color-magnitude 
diagram of this galaxy. The lower Figure shows the predictions of a model with two 
separate bursts which seems to fit better the data. 
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Fig. 46. Observed and predicted [Mg/Fe] vs. [Fe/H] relation for the dSph Draco, 
Ursa Minor and Sextans. The different curves refer to different SF efficiencies (csf) 
expressed in Gyr -1 , which are equivalent to the quantity v. Figure from Ikuta & 
Arimoto (2002). 



• a strong galactic wind develops when the thermal energy of the gas equates 
the binding energy of the gas. The rate of gas loss is assumed to be several 
times the SFR, as in eq. (27), with typical values of A (wind parameter) 
between 5 and 15. 

The condition for the onset of the wind is written as: 

(Eth)lSM > Esgas (37) 

namely that the thermal energy of the gas is larger or equal to its binding 
energy. The thermal energy of gas due to SN and stellar wind heating is: 

(Eth)iSM = E t h SN + Eth w (38) 
with the contribution of SNe being: 

EthsN — / e SNRsN{t )dt , (39) 
Jo 
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while the contribution of stellar winds is: 

E thw = <p{m)ip(t)e w dmdt (40) 

JO J12 

with €sn — Vsn^o an d e o = 10 51 erg (typical SN energy), and e w = tj w E w 
with E w = 10 49 erg (typical energy injected by the stellar wind of a 20M Q 
star taken as representative). rj w and tjsn & re two free parameters and 
indicate the efficiency of energy transfer from stellar winds and SNe into 
the ISM, respectively, quantities still largely unknown. It is assumed that 
e w — 0.03 for the stellar winds, and that €sn = 0.03 for Type II SNe and 
csn = 1-0 for Type la SNe, as suggested by Recchi et al. (2001). The total 
mass of the galaxy is expressed as M to t{t) = M*(t) + M gas (t) + Mdark{t) 
with M^t) — M*(t) + M gas (t) and the binding energy of gas is: 

E Bgas (t)=W L (t) + W LD (t) 

with: 

w L (t) = _ . 5g jWW) 

which is the potential well due to the luminous matter and with 

„ T ... „ M gas {t)M dark 

W L D(t) = —Gwld — — (43) 

which represents the potential well due to the interaction between dark 
and luminous matter, where wld ~ + 1.375), with S = r^/ro, 

being the ratio between the galaxy effective radius and the radius of the 
dark matter core. The typical model for a dSph starts with an initial 
baryonic mass of 1O 8 M0 and ends up, after the wind, with a luminous 
mass of ~ 1O 7 M . The dark matter halo is assumed to be ten times larger 
than the luminous mass but diffuse (S = 0.1). The galactic wind in these 
galaxies develops after several Gyr from the start of SF, according to the 
different assumed SF efficiency. Soon after, but not immediately, the wind 
has started, the SF decreases very strongly until it halts completely. 

• The IMF is that of Salpeter (1955) for all galaxies. 

• Each galaxy is supposed to have formed by infall of gas clouds of primordial 
chemical composition, on a timescale not longer than 0.5 Gyr. 

In Figure 47 we show the [a/Fe] ratios for different a-elements and for 
different efficiencies of SF, as predicted by the standard model of Lanfranchi 
& Matteucci (2003). 

As one can see, the [a/Fe] ratios show a clear change in slope followed by a 
steep decline, in agreement with the data. The change in slope corresponds to 
the occurrence of the galactic wind which starts emptying the galaxy of gas. In 
such a situation the SF starts to decrease as does therefore the production of 
the a-elements from massive stars, whereas Fe continues to be produced since 



(41) 
(42) 



Chemical evolution of the Milky Way and its Satellites 65 




[Fe/H] 

Fig. 47. Observed and predicted [a/Fe] versus [Fe/H]. The different lines refer to the 
"standard model" with different SF efficiency v going from 1 (dashed-dotted lines) 
to 0.01 Gyr^ 1 (continuous lines). The points represent stars in different dSphs: 
Sagittarius (open triangles), Draco (filled hexagons), Carina (filled circles), Ursa 
Minor (open hexagons), Sculptor (open circles), Sextans (filled triangles), Leo I 
(open squares) and Fornax (filled squares). Figure and references from Lanfranchi 
& Matteucci (2003). 

its progenitors have long lifetimes. This produces the steep slope: the low SF 
efficiency and the wind, which decreases furtherly the SF. In this situation, the 
time-delay model predicts an earlier and steeper decline of the [a/Fe] ratios, 
as we have already discussed. 

In LM04, the histories of star formation of specific galaxie, as suggested 
by their color-magnitude diagrams, were taken into account, and they devel- 
oped models for six dSphs: Carina, Ursa Minor, Sculptor, Draco, Sextans and 
Sagittarius. In Table 1 we show the assumed SF formation histories and model 
parameters. In particular, in column 1 are the galaxy names, in column 2 are 
the initial baryonic masses, in column 3 are the SF efficiencies, in column 4 are 
the wind parameters, in column 5 are the numbers of SF episodes, in column 
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6 are the times at which the SF episodes start, in column 7 are the durations, 
in Gyr, of the SF episodes and in column 8 the assumed IMF. 

Table 1. Models for dSph galaxies. M™l tlal ; s the baryonic initial mass of the 
galaxy, v is the star-formation efficiency, A is the wind efficiency, and n, t and d are 
the number, time of occurrence and duration of the SF episodes, respectively. 
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In Figures 48, 49, 50, 51 and 52we show the predictions for specific dSphs 
by LM04. As one can see, the [a/Fe] data in the dSphs are well reproduced 
and in particular the steep decline of the [a/Fe] ratio is well reproduced. 
This steep decline is due again to the low efficiency SFR, a feature common 
also to the other models, coupled with a strong and continuous galactic wind 
which gradually empties the galaxies of gas. In the previous models either 
the galactic wind was not present or it was assumed instantaneous or not as 
strong as in LM04, thus predicting a flatter slope for the descent of the [a/Fe] 
vs. [Fe/H], as we have already seen. 

Lanfranchi et al. (2006) computed also the expected abundances of s- and 
r- process elements in dSphs, by adopting the same nucleosynthesis prescrip- 
tions used for the chemical evolution of the Milky Way. In particular, they 
adopted the prescriptions of Cescutti et al. (2006) for Ba, Y, La, Sr and Eu: 
Ba, Sr , La and Y are mainly s-process elements produced on long timescales 
by low mass stars (1 — 3M Q ), but they have also a small r-process component 
originating in stars in the mass range 12-30M Q . The Eu instead is considered 
as a pure r-process element produced in the stellar mass range 12-30M Q . 

In Figures 53, 54, 55 and 56 we show the predictions of Lanfranchi et al. 
(2006) for s-and r-process elements in dSphs compared with the available data. 
Also in this case the agreement looks good, although more data are necessary 
before drawing firm conclusions. The general tendency for the a-elements in 
dSphs is to be less overabundant relative to Fe and the Sun than the stars of 
the solar vicinity with the same [Fe/H]. This is due to the lower SFR in dSphs 
(the effect is increased by the galactic wind) which acts to shift the curve 
[a/Fe] vs. [Fe/H] for the solar vicinity towards left in the diagram, whereas 
a stronger SF than in the solar neighbourhood moves the solar vicinity curve 
toward right in the diagram (see Figure 33). The same holds for s- and r- 
process elements: in this case, since [s/Fe] vs. [Fe/H] first increases sharply 
at low metallicities and then it flattens at higher ones (the opposite of what 
happens for the a-elements), the dSphs show a higher [s/Fe] than the stars 
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Fig. 48. Observed and predicted [a/Fe] vs. [Fe/H] relation for the galaxy Carina. 
The different lines represent models with different SF efficiency. The continuous line 
represents the best model and corresponds to the efficiency v — O.lGyr^ 1 . Figure 
from LM04. 

in the solar vicinity at the same [Fe/H]. This shows again the effect of the 
time-delay model. 

The differences between the [s/Fe] vs. [Fe/H] in dSphs and the Milky Way 
are shown in Figures 57, where data and predictions of Ba and Eu for Sculptor 
are compared with data and predictions of Ba and Eu in the solar vicinity. 

Finally, another important constraint for models of galactic chemical evo- 
lution is represented by the stellar metallicity distribution. In Figures 58 we 
show the predictions for the stellar metallicity distribution of Carina com- 
pared with the observed one and the agreement is very good. The observed 
distribution is from Koch et al. (2006) who measured the metallicity of 437 
giants in Carina by means of Ca triplet and then transformed it into [Fe/H] 
through a suitable calibration. In Figures 58 we also show the comparison 
between the stellar metallicity distribution in Carina and the G-dwarf metal- 
licity distribution in the solar vicinity. As one can see, the Carina distribution 
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Fig. 49. Observed and predicted [a/Fe] vs. [Fe/H] relation for the galaxy Sculp- 
tor. The different lines represent models with different SF efficiencies (y — 
0.05, 0.2, 0.5Gyr _1 ). The continuous line represents the best model {y = 0.2G?/r _1 ). 
Figure from LM04. 



lies in a range of smaller metallicities due to the lower efficiency of SF assumed 
for this galaxy. 

To this purpose, owing a word of caution is appropriate: in fact, the Ca 
triplet in principle traces the abundance of Ca and not that of Fe and we 
know that Ca and Fe evolve in a different way since Ca is mainly produced in 
Type II SNe, whereas Fe is produced mainly in Type la SNe. This different 
evolution of Ca and Fe leads, in the Koch paper, to an uncertainty of 0.2 dex. 
Besides that, the globular clusters which serve as calibrators for obtaining 
[Fe/H] lie in the range -2.0 - -1.0 dex, whereas Koch's data extend down to 
lower metallicities. 

The good fit of the stellar metallicity distribution indicates that both the 
assumed history of SF and the IMF are close to reality. 

LM04 predicted the stellar metallicity distribution for all the six dSphs 
and while for Carina the agreement is quite good for Sculptor they cannot 
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Fig. 50. Observed and predicted [a/Fe] vs. [Fe/H] relation for the galaxy Ursa 
Minor. The different lines represent models with different SF efficiencies. The con- 
tinuous line represents the best model (y = 0.2Gj/r _1 ). Figure from LM04. 

reproduce the bimodal stellar distribution recently suggested by Tolstoy et al. 
(2004) and shown in Figure 59, since their model is a one-zone model. In Figure 
60 are shown the predictions of LM04 for the Sculptor galaxy: it is clear from 
this Figure that to reproduce the two different stellar populations, one has to 
assume a multizone model possibly with different efficiencies of SF. Kawata 
et al. (2006) explained the bimodality of stellar populations in Sculptor as a 
consequence of dissipative collapse which produces higher metallicities at the 
center of the galaxy. 

7 What have we learned about dSphs? 

From the study of the chemical evolution of dSphs and the Milky Way we can 
derive the following conclusions: 
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Fig. 51. Observed and predicted [a/Fe] vs. [Fe/H] relation for the galaxy Draco. 
The different lines represent models with different SF efficiencies. The continuous 
line represents the best model (y — O.OSGyr^ 1 ). Figure from LM04. 

• By comparing the [a/Fe] ratios in the Milky Way and in dSphs of the 
Local Group we can conclude that these systems had different histories of 
SF. 

• The [a/Fe] ratios in dSphs are always lower than in the Milky Way at the 
same [Fe/H]. This is a consequence of the time-delay model which predicts 
this behaviour for systems which suffered a lower star formation than the 
solar vicinity. 

• The occurrence of strong galactic winds or gas loss in general is necessary 
to keep the SF low and it produces the steep decrease of the [a/Fe] ratio 
observed in dSphs (see Lanfranchi & Matteucci 2007). 

• Good agreement is found both for [s/Fe] and [r/Fe] abundance ratios. The 
[s/Fe] ratios are predicted to be higher than the same ratios in Milky Way 
stars with the same [Fe/H]. This is again a consequence of the time-delay 
model. 
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Fig. 52. Observed and predicted [a/Fe] vs. [Fe/H] relation for the galaxy Sagit- 
tarius. The different lines represent models with different SF efficiencies. The con- 
tinuous line represents the best model(^ = S.OGyr^ 1 ) . Here the efficiency of SF is 
relatively high but the strong galactic wind makes it effectively much lower. Figure 
from LM04. 



• The dSphs of the Local Group contain very old stars but they suffered 
extended periods of SF, far beyond the reionization epoch. This is sug- 
gested both from the color-magnitude diagrams of these galaxies and from 
the level of the abundances of s-process elements such as Ba, which could 
not have been observed if the SF had stopped at the reionization epoch 
(Fenner et al. 2006). 

• All the previous conclusions suggest that it is unlikely that the dSphs have 
been the building blocks of the Milky Way as predicted by current CDM 
models (see review by Geisler et al. 2007 for a detailed discussion of this 
point). Robertson et al. (2005) studied the formation of the Galactic stellar 
halo by means of different accretion histories for the dark matter halo of 
the Milky Way in the framework of the ACDM model. They concluded, on 
the basis of the [a/Fe] ratios in Galactic halo stars, in dwarf irregulars and 
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Fig. 53. Predicted and observed [s,r/Fe] vs. [Fe/H] for the galaxy Carina. The 
model is from Lanfranchi, Matteucci & Cescutti (2006a), where the references for 
the data can be found. 



dSphs, that it is more likely that the Galactic dark matter halo was formed 
by an early accretion of dwarf irregular galaxies, which formed stars for a 
short time and then were destroyed. Concerning dSphs, they suggest that 
their chemical abundances should have been affected by galactic winds and 
that they should have been accreted and destroyed over the entire Milky 
Way lifetime. 



8 Other spirals 

In this section we will briefly describe the properties of other spirals of the 
Local Group for which we have enough chemical information. 
In external spirals we can measure: 

• the SFR mainly from H a emission and this suggests a correlation between 
the SFR and the total surface gas density, as discussed previously. 
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Fig. 54. Predicted and observed [s,r/Fe] vs. [Fe/H] for the galaxy Sculptor. The 
model is from Lanfranchi, Matteucci & Cescutti (2006a), where the references for 
the the data can be found. 



Abundance gradients are also found in disks of local spirals (see Garnett 
et al. 1997): in particular, it is found that abundance gradients, expressed 
in dex/Kpc, are steeper in smaller disks but the correlation disappears 
if they are expressed in dex/ Rd (where Rd is the disk scalelength) . It is 
interesting to note that a universal gradient slope per unit scale length may 
be explained by viscous disk models (e.g. Lin & Pringle, 1987; Sommer- 
Larsen & Yoshii, 1989). Further information on gradients is that they are 
flatter in disks with bars: probably the bar induces radial flows which can 
wash-out the abundance gradients if their velocity is high enough (sec 
Tinsley, 1980). 

The gas distribution: one finds differences in the gas distribution along the 
disk betwen field and cluster galaxies, these latter being subject to ram 
pressure stripping. 

Integrated colors of galaxy disks give information on the distribution of 
stellar populations along the disks. Several authors (Josey & Arimoto, 
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Fig. 55. Predicted and observed [Ba/Eu] vs. [Fe/H] for the galaxy Sculptor. The 
model is from Lanfranchi, Matteucci & Cescutti (2006a), the data and the figure 
are from Geisler et al. (2007). 



1992; Jimenez et al. 1998; Prantzos & Boissier, 2000) have suggested that 
color gradients, as well as metallicity gradients, can be reproduced by 
assuming an inside-out formation for disks, as has been suggested for the 
disk of the Milky Way. 
• In Figure 61 we show the results of a paper by Boissier et al. (2001). 
They conclude that more massive disks are redder, more metal rich and 
more gas-poor than smaller ones. On the other hand, their estimated SF 
efficiency (defined as the SFR per unit mass of gas) seems to be similar 
among different spirals: this leads them to conclude that more massive 
disks arc older than less massive ones. The various quantities in Figure 61 
are plotted as functions of the disk circular velocity which is a measure of 
the dark matter halo of each galaxy. The various curves are obtained by 
varying the spin parameter A, expressed as a X/Xmw, where Xmw refers 
to the Milky Way. In fact, in the framework of semi- analytical models of 
galaxy formation the evolution of galactic disks can be described by means 
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Fig. 56. Predicted and observed [s,r/Fe] vs. [Fe/H] for the galaxy Ursa Minor. The 
model is from Lanfranchi, Matteucci & Cescutti (2006a), where the references for 
the the data can be found. 



of scaling laws calibrated on the Milky Way with V c and A as parameters 
(e.g. Mo et al. 1998). 

8.1 Chemical models for external spirals 

Several models of chemical evolution of Local Group spirals have been devel- 
oped in the past years (e.g. Diaz & Tosi, 1984; Molla et al. 1996; Chiappini 
et al. 2003). Diaz & Tosi (1984) first modeled the chemical evolution of M31, 
M33, M83 and M101. Molla et al. (1996) modeled several spirals of the Local 
Group (M31, NGC300, M33, NGC628, NGC3198, NGC6946). In Figure 62 
we show the predictions, compared to observations, of the Molla et al. model 
for M31. 

As an another example of abundance gradients and gas distribution in 
a local spiral galaxy we show in Figure 63 the observed and predicted gas 
distribution and abundance gradients for the disk of M101. In this case the gas 
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Fig. 57. Comparison between data and predictions for Ba and Eu in Sculptor and 
in the Milky Way. Model for the Milky Way: continuous lines. Model for Sculptor: 
dashed lines. Data for the Milky Way: triangles. Data for Sculptor: full squares. 
Models from Lanfranchi & al. (2007) where the references to the data can be found. 





Fig. 58. Upper panel: stellar metallicity distribution for Carina. Model from Lan- 
franchi et al. (2006b). The assumed SF efficiency is v — 0.15G?/r _1 and the wind 
efficiency is A = 5. Two different histories of SF have been tested here: the one of 
Dolphin (2002) (continuous line) and that of Rizzi et al. (2003) (long dashed line), 
but this does not produce important differences in the results. The main difference 
between the two histories of SF is the number of bursts (3 in Dolphin and 4 in Rizzi 
et al.). Data from Koch et al. (2006). Lower panel: predicted stellar metallicity dis- 
tribution for Carina compared with the predicted G-dwarf metallicity distribution 
in the solar neighbourhood (dashed line) . 
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Fig. 59. Observed stellar metallicity distribution for Sculptor. Data and figure are 
from Tolstoy et al. (2004): all the stars of Sculptor are indicated by the dotted line. 
The central stars are those indicated by the lower histogram with continuous line, 
whereas the stars beyond R > 0.2 Kpc are indicated by the upper histogram with 
continuous line. 



distribution and the abundance gradients are reproduced with systematically 
smaller timescales for the disk formation relative to the MW (M101 formed 
faster), and the difference between the timescales of formation of the internal 
and external regions is smaller (jMim = 0.75r(Kpc) — 0.5 Gyr, Chiappini et 
al. 2003) compared to eq. (34). 

Therefore, in the Chiappini et al. (2003) paper it is suggested that the 
fact that more luminous spirals ( e.g. M101) tend to have shallower abun- 
dance gradients than less luminous ones can be interpreted as due to a faster 
formation (down-sizing) of large spirals, as also hinted at by the results of 
Boissier et al. (2001). 

In summary, from the available studies of spirals in the Local Group we 
can suggest the following: 
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Fig. 60. Observed and predicted stellar metallicity distribution for Sculptor. The 
data (the dotted line of Figure 59) are represented by the histogram (long dashed). 
The models are from LM04: the solid line represents the best model, whereas the 
dashed lines represent models with higher (the curve on the right) or lower (the 
curve on the left) SF efficiency. 



• The disks of spirals have all formed inside-out and the more massive disks 
have formed faster than the less massive one. 

• This translates into a faster gas accretion rate and consequently into a 
faster SFR. 

• In other words, the most massive disks are also the oldest, a conclusion 
which is not in line with classical predictions from CDM models. 



9 Cosmic chemical evolution 

With the name cosmic chemical evolution we indicate the chemical evolution 
taking place in comoving volumes large enough to be representative of the 
whole universe (Pei & Fall 1995). The evolution can be described in terms of 
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Fig. 61. Distribution of the (B-H) color, the 12+log(0/H), the gas fraction and the 
SFR efficiency (the SFR per unit mass of gas) versus the circular velocity of disks 
of local spirals. The different lines refer to different values of the spin parameter A. 
Figure from Boissier et al. (2001) where the references to the data can be found. 
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Fig. 62. Predicted and observed abundance gradients in M31. Figure and models 
are from Molla et al. (1996). 
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Fig. 63. Upper panel: predicted and observed gas distribution along the disk of 
M101. The observed HI, H2 and total gas are indicated in the Figure. The large 
open circles indicate the models: in particular, the open circles connected by a 
continuous line refer to a model with central surface mass density of W00Mqpc~ 2 , 
while the dotted line refers to a model with 800Mopc~ 2 and the dashed line to a 
model with 600Mqjjc - . Lower panel: predicted and observed abundance gradients 
of C,N,0 elements along the disk of M101. The models are the lines and differ for a 
different threshold density for SF, being larger in the dashed model. All the models 
are by Chiappini et al. (2003). 
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comoving densities of gas and stars f2 gas and f2 s tars, both measured in units 

of the present critical density (p c — §^j) and the mean abundance of heavy 
elements in the ISM, Z, including dust. Under IRA we can write, following 
Pei & Fall: 

^gas ~t~ ^stars — ^ f (44) 

and 

HgasZ - VZ^ stars = (Zf - Z)Q f (45) 

where the dots represent differentiation with respect to the cosmic time. 
The term Qf can represent either the infall or the outflow rate according to 
its sign. For the closed-box model J?/ = and its solution is: 

Z = -yzln{Qgas/ Mgas^ ) (46) 

with flgas^, being the gas comoving density at some suitably high redshift 
when there are still no stars and heavy elements. 

By means of these equations Pei & Fall followed the evolution of DLAs 
(the quasar absorbers). They expressed the quantity f} gas in terms of the 
observable properties of DLAs. They assumed IRA. 

Since all these cosmic quantities refer to an unitary volume of the universe 
which contains galaxies of all morphological types, 

Calura & Matteucci (2004) proposed another approach to the cosmic chem- 
ical evolution, which takes into account galaxies of different morphological 
type. They computed the cosmic chemical enrichment of the universe by means 
of detailed models of chemical evolution of galaxies of all morphological types, 
relaxing IRA and assuming for each galaxy type a different history of SF, as 
discussed in the previous sections. They defined the comoving cosmic density 
of stars and gas for galaxies of different morphological type (ellipticals, spirals 
and irregulars) as: 

P*,k = PB.k ■ (M*/L) B ,k (47) 

for the stars and 



Pg.k = PB.k ■ {M g /L) B>k (48) 

for the gas. The quantities (M*/L)_e,fe and (M g /L)B,k are the predicted M/L 
ratios for stars and gas, respectively. Ls,k is the blue luminosity for each 
galaxy type (k indicates the morphological type) and ps,k is the comoving 
luminosity for a given galaxian morphological type. 

They computed the mean mass weighted metallicity of galaxies by sum- 
ming the metallicities predicted for the different morphological types as: 



X)fc Pg,k Zg.k + X)fe P*,k Z* v k 



, 7 . _ Z^fc h'gM ^gM T /_, k y*M ^*,k ,.„n 

*^ ^galaxies ^ ^ f ^ , ^ \^ ) 
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and obtained: 



< Z galaxtes >= 0.0175 = 0.9Z Q (50) 

where Z Q = 0.02 and with 56%, 42% and 2% of metals produced in el- 
lipticals, spirals and irregulars, respectively. Therefore, the conclusion is that 
the average metallicity in galaxies is almost solar and that most of the metals 
in the universe have been produced by elliptical galaxies. They also predicted 
the average [O/Fe] ratios for each galaxy type both in the gas and in stars. 

In particular: 

[0/Fe]*,BjM pi = 0.4 dcx, [0/Fe] gas , E Uip = -0.33 dex, [O/Fe]*, spiral = 0.1 
dex and [0/Fe] 9OS) sp iro j = O.Oldcx. 

Then, they computed the metallicity in the intergalactic medium (IGM) 
by considering all the metals ejected by galaxies (mainly ellipticals) into the 
IGM: 

. ry . &ZJGM 

< Z IGM >= -p. (51) 

Mb, IGM 

with fib, igm — ^b — ^b.* — ^b. gas = 0.0753 being the baryonic density of 
the IGM and fl h = 0.02/i~ 2 (from WMAP, Spergel et al. 2003, 2007) being 
the total baryonic content of the universe. Therefore, they obtained: 

< Z IGM >= 6.54 x 10~ 4 = O.O3Z . (52) 

Finally, they computed the average metallicity of the universe by account- 
ing for all the metals produced in galaxies over the lifetime of the universe: 

< Z universe > = E*L%i = = 0.0017 = O.O9Z (53) 

J4 lit, 

where flz,tot represents the sum of all the metals produced in all galaxies and 
fib represents the the total amount of baryons in the universe. 

In summary, the mean metallicity inside galaxies of all morphological types 
is almost solar, whereas the mean metallicity of the universe is roughly 1/10 
solar. 

Calura & Matteucci (2003) computed also cosmic Type la SN rates 
(SNR cosm ), expressed in SNu (number of SNe/10 10 LB Q per century). In 
particular, they took into account the contribution of all galaxy types in the 
following way: 

SNR cosm (z) - E t?T fc(Z) ( 54 ) 

where SNR cosm can represent the Type II, Ib/c, la SNe. The sums are 
over all the galactic morphological types, Ls k is the total blue luminosity of 
the kth morphological type. In order to compute the SNR cosm for each galaxy 
they assumed a SN model progenitor and a cosmic SFR, calculated as: 
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P* = ^2pbM z ) ■ ( M */ L )B,k(z) -i>(z)k, (55) 
k 

where pB,k{z) and {M*/ L)b,u{z) have been already defined and ip(z)k repre- 
sents the history of SF of a galaxy of kth morphological type. They assumed 
the SF histories of Figure 41. In Figures 64 and 65 we show some examples 
of predicted cosmic SN rates by adopting the same cosmic SFR (eq. 55) but 
different assumptions about the Type la progenitor model. 

As one can see from these Figures, the best DTD appears to be the one of 
Mannucci et al. (2006), although this conclusion is based only on the fit of the 
rates in radio galaxies. Concerning the cosmic Type la SN rate (panel b) the 
agreement is good for both DTDs except for the point at the highest redshift 
which is highly uncertain. We need more detections of SNe la at high redshift 
before drawing any conclusion on the high redshift cosmic Type la SN rate. 
In conclusion, the study of the Type la SN cosmic rates is very important to 
impose constraints on the SN progenitor model, the histories of SF in galaxies 
and, last but not least, the Hubble diagram. 
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Fig. 64. Theoretical cosmic Type la SN rates compared with observational data 
from Mannucci et al. (2006). The progenitor model adopted assumes the delay-time 
distribution suggested by Mannucci et al. (2005;2006) (panel a), whereas the cosmic 
SFR is the one of Calura & Matteucci (2003) (panel b). In panels c) an d) are shown 
predictions and data for the Type la SN rate per unit galactic mass (SNuM) versus 
color and radio flux in radio galaxies, respectively. For the references about the data 
see Mannucci et al. (2006). Figure adapted from Mannucci et al. (2006). 
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Fig. 65. Theoretical cosmic Type la SN rates compared with observational data 
from Mannucci et al. (2005;2006). The progenitor model adopted assumes the delay- 
time distribution suggested by Matteucci & Recchi (2001) (panel a), whereas the 
cosmic SFR is the one of Calura & Matteucci (2003) (panel b). In panels c) an d) are 
shown predictions and data for the Type la SN rate per unit galactic mass (SNuM) 
versus color and radio flux in radio galaxies, respectively. For the references about 
the data see Mannucci et al. (2006). Figure adapted from Mannucci et al. (2006). 
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